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Chapter 1 

Introduction 



Quantum impurity models describe interactions between some local degrees of freedom (e.g. a 
spin) and a continuum of non-interacting fermionic or bosonic states. The investigation of quan- 
tum impurity models is a starting point towards the understanding of more complex strongly 
correlated systems, but quantum impurity models also provide the description of various cor- 
related mesoscopic structures, biological and chemical processes, atomic physics and describe 
phenomena such as dissipation or dephasing. Prototypes of these models are the Anderson 
impurity model, or the single- and multi-channel Kondo models. The first two models are clas- 
sic examples of Fermi liquid models, while the multi-channel Kondo model is the most basic 
example of a non-Fermi liquid system, and as such, it serves possibly as the simplest realization 
of a quantum critical state. 

The solution of these models for low energies was a major issue in theoretical condensed matter 
research and led to the development of various non-perturbative techniques. [The interested 
readers are referred to the seminal book of Hewson [T] and to the extensive review of Cox and 
Zawadowski [2J.] However, despite the extensive effort, many of the methods developed are 
uncontrolled, while others can be applied only to a subclass of models or to restricted regions 
of the parameter space. Wilson's numerical renormalization method, originally developed for 
the Kondo model, remained possibly the most reliable method to study dynamical correlations 
of generic quantum impurity models as well as their thermodynamic properties and finite size 
spectra. It is still one of the most popular methods to study quantum impurity models. 
Although Wilson's NRG has been used in its original form for a longtime, a number of new de- 
velopments took place recently: First, a spectral sum-conserving density matrix NRG approach 
(DM- NRG) has been developed [21 0] , which has very recently been generalized for non-Abclian 
symmetries [5]. We remark that using symmetries as much as possible is necessary in many 
cases to perform accurate enough calculations. NRG has been also restructured as a matrix 
product state approach [5J. These new developments not only made NRG much more reliable 
than Wilson's original method [5J [7] , but they made it possible to extend NRG to study non- 
equilibrium phenomena [HI E] , and opened the way to use methods familiar from the density 
matrix renormalization group (DMRG) approach [TUj . 

In this manual we do not intend to give a complete description of the NRG machinery. Rather, 
we introduce some of the basic concepts that are needed to use NRG and to use the code we 
provide for quantum impurity problems of interest. If you want to understand in detail, how 
NRG and DM-NRG work, we advise you to consult Wilson's original work [TTj, [T2] , and the 
references listed above. 

The code we describe in this manual is a free density matrix numerical renormalization group 
(DM-NRG) code, which can be downloaded from the site http : //www.phy .bme .hu/^dmnrg. 
This code is a flexible NRG code, which uses user-defined non-Abelian symmetries dynamically, 
computes spectral functions, expectation values of local operators for user-defined impurity 
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models. The code can use a uniform density of states as well as a user-defined density of states. 
The current version of the code assumes fermionic bath's. It uses any number of U(l), SU(2) 
charge SU (2) or Z% symmetries, but the interested user can teach the code other symmetries 
too. The code runs using a simple input file. We provide several example input files with the 
code as well as a few Mathematica files with which these input files can easily be constructed. 
An energy spectrum analyzer utility is also provided with the code to study finite size spectra 
too. 
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Chapter 2 



A short introduction to 
DM-NRG and to the use of 
symmetries 



In this chapter, we give a short overview of Wilson's numerical renormalization group (NRG) 
and the density matrix NRG (DM-NRG). As already mentioned in the introduction, here we 
introduce only the basic concepts that are needed to use NRG and to use the code we provide 
for quantum impurity problems, but we do not attempt/intend to give a complete overview of 
the existing literature. To learn more details about NRG and DM-NRG, we recommend to read 
Wilson's original work [TTJ [T2] . 

2.1 Wilson's NRG 

2.1.1 The simplest example: The Kondo model 

Before introducing the general concepts, let us discuss the simplest possible example, the Kondo 
model. The Kondo model consists of a spin S interacting locally with a non-interacting con- 
duction electron sea, 

H-Kondo = — ^0> iO .'VV' + #cond • (2-1) 
(j, a' 

Here J denotes the Kondo coupling, ip^ creates a conduction electron at the impurity site, and 
a a ,a-' is the vector of Pauli spin matrices. The term, -ff C ond, describes the conduction electron 
bath, and its specific form is not very important for us. 

Wilson's mapping and iterative diagonalization 

From the point of view of local dynamics, H con d contains a lot of redundant information: In fact, 
since our fermions are non-interacting, the local density of states q(u>), i.e. the spectral function 
of the unperturbed Green's function, (w), is generated by -ffcond with J = 0. Note that 

g(u>) determines completely the correlation functions of 4>t & n d the spin dynamics even for the 
interacting system, J ^ 0. The ingenious idea of Wilson was to discretize logarithmically g(uj) 
using a discretization parameter A > 1, and thus map approximately the Hamiltonian (|2.1| to 
a semi- infinite chain (see Fig. 12.11 and Fig. 12. 2[) : 
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Figure 2.1: Discretization of the conduction band density of states on a logarithmic mesh and 
mapping onto the Wilson chain (see also Fig. I2.2|) . 



h%:Lt = ^E fl*3°0>h*> + E E^« ti,*f»,° + E E*» + ■ ( 2 - 2 ) 

<7,a"' n— cr n— <r 

Here the spin interacts only with the fermion /q a at the end of the Wilson chain. The operator 
/q a creates a conduction electron at the impurity site, and is essentially identical to the operator 
ip\. The on-site energies and the hoppings t n depend solely on g(u>) and can be determined 
recursively [11] . Our code takes care of this part of the work: it determines numerically these 
constants if a density of states g(cu) is provided (see section 14.4.31 and the description of the 
utility he). 

For a flat and symmetrical density of states, g(ui) = l/(2D), one has £ n = and the t n 's can 
be determined analytically ( t„ ~ hr n / 2 ) [TTj. Longer and longer chains give more and more 
accurate description of the infinite chain. This observation led Wilson and his co-workers to 
solve the Hamiltonian (|2.2[) iteratively. Introducing the operator 

H„ = ^^E fo,<T°<r,<r' fo,<r' + E E fm,afm,a + E E tm (fL,afm+l,a + h.C.) . (2.3) 
cr, cr' m— cr rn—0 a 

one has the obvious recursion relation, 

Hn+i = H n + T nJl+ i + Hn+l , (2.4) 

with the notation, T n , n +i = *nX)a(/i,o-/n+i,<r + h.c.) and H„+i = Z) CT fn+i fl+i.afn+i,* (see 
also Fig. l2.2p . Then Wilson's procedure consists of constructing from the low-energy eigenstates, 
\u) n , of the operator H n approximate eigenstates, |u) rl +i, of the operator H n+ i. To do this, 
one takes a definite number of the states \u) n with the lowest energies and generates new states 
from them by first adding an empty site, \u) n — > |u) n +i and then using the operators „ to 
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Figure 2.2: Wilson chain of length N. The impurity site is at the origin and is represented by 
a square, while the sites are represented as circles. For the impurity the Hamiltonian 7io is 
identical with Hq. As we add sites the iterative construction of the Hamiltonian at a given site 
H n is pictured. 



create new states as 



/ \u)n+l \ 

f n +l,l\ u )n+l 
\ fn+l,tfn+l,l' U ' n+1 / 



(2.5) 



Then one diagonalizes the Hamiltonian H n+ i in this new basis set. To construct the matrix 
elements of H n+ i in this new basis, one needs the following information 

1 . The eigenvalues E™ of H n , 



2. The matrix elements of 7i n +i between the local states, \fi), constructed from the vacuum 
state |0) as 

/ 1°) \ 

4 + i,tl°>. 

V /l + i, T 4 + i,JO) J 



(2.6) 



3. The matrix elements of a between the local states, \fi), and 

4. The matrix elements of f\ a between the block states \u) n . 

Diagonalizing the Hamiltonian H n+ i one then obtains the new eigenstates |m)„+i, their eigen- 
values, E? +1 , and one can trivially compute the matrix elements of a between these states 
too, to proceed to the next iteration. 

If one is to compute the spectral function of a local operator A acting at site n = 0, then one 
has to keep track of the matrix elements „(u|A|u) n of this operator, too. Already the finite size 
spectrum, i.e. the spectrum of H n contains a lot of precious information However, once 
the matrix elements and the approximate spectrum of the Hamiltonian is at hand, one can go 
ahead and also compute thermal expectation values or spectral functions from it [T3J [14] , and 
thus determine completely the properties of a quantum impurity problem. 
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Symmetries in the Kondo model 



The total spin operator 

St = S+-J2 E A^°,°'fn,*' (2-7) 

n=0 (T,<7'e{T,i} 

as well as the charge operator, 

Q = E E Arfn,* - 1 ( 2 - 8 ) 

n=0 \<r={T,!} / 

commute with Eq. (|2.2[) . This implies that the eigenstates of the Hamiltonian can be classified 
as multiplets. Every multiplet u will be characterized by a charge quantum number, Q{u) and 
a spin quantum number, S(u). We shall refer to these as representation indices. Internal states 
within a multiplet are labeled by the z-component of the spin, that we shall refer to as labels. 
Similarly, operators can also be characterized by quantum numbers. To give an example, in this 
simple case, the components of {f n p f n ,i} and {/^ , , — /n,f} transform under spin rotations as 
the components | ± 1/2) of a spin Sf = 1/2 state, and have charges Qf, =1/2 and Qf = —1/2, 
respectively. The three operators formed from the impurity spin, S m = {—S + /\^2, S z , S~ /V2} 
transform, on the other hand, as the three components of a spin Ss = 1 state, \m = 0, ±), while 
they have trivially charge Qs = 0. The operators discussed before provide simple examples 
of irreducible tensor operator multiplets {A}. Notice the non-trivial prefactors in the previous 
examples. These prefactors must always be worked out carefully, since the smallest sign mistake 
may alter the NRG results significantly. 

In the presence of symmetries, a powerful theorem, the Wigner-Eckart theorem tells us that, 
in order to compute the matrix elements of a member of an operator multiplet {A} between 
two states within multiplets u and v a single matrix element is needed, the so-called reduced 
matrix element [5] : 

(u\\A\\v) . (2.9) 

Wilson's procedure thus gets a bit modified. First of all, one needs to classify not only the 
block states (multiplets) but also the local states fi added at site n + 1 by symmetries. Along 
the iteration, one constructs from these new states using group-theoretical methods (Clcbsch- 
Gordan coefficients), which are already eigenstates of Sj, and S z . Apart from this twist, the 
discussion of the previous subsection is still valid, and gets only slightly modified. We can thus 
make the following statement: To construct the (reduced) matrix elements of (i||i? n+ i||j) in 
this basis, one needs the following information 

1. The eigenvalues E™ of H n , 

2. The matrix elements (/z||W n +i||i / ) between the local multiplets. 

3. The reduced matrix elements, (mII/I+iH^), 

4. The reduced matrix elements n (u\ \f\ \ \v) n between the block multiplets, \u) n . 

5. Finally, to compute a spectral function of an operator multiplet A, one also needs „ (u\ \ A\ \v) 

At this point, the reader does not need to know in detail, how the construction of H n+ \ takes 
place, since the code does that automatically, but for more details, the interested reader is 
referred to Ref. [S]. 
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Table 2.1: Block states (left) and local/added states (right) for the single channel Kondo model 
when U c harge{^) x U sp i n (l) symmetry is used. The first block states are formed from the 
impurity spin and the conduction electron at site n = of the Wilson chain. 



Kondo model with U c harge(l) X U sp i n (l) symmetry 



Let us first illustrate the concepts introduced above through the simplest case, when every 
symmetry used is Abelian. In this case every multiplet consists of a single state, which is 
characterized by some quantum numbers. Let us thus consider the Kondo model with two 
£7(1) symmetries, generated by the ^-component of the total spin, S z and the charge operators, 
Q. Now S z and Q, provide the quantum numbers according to which the multiplcts of the 
Hamiltonian arc classified. The classification of the states in this case is rather trivial and 
the computation of the reduced matrix elements is also simple, since all the Clebsch-Gordan 
coefficients are 1 or 0. However, without spin SU sp i n (2) symmetry, the hoppings of ^ and 

^ are not related by symmetry. There are therefore two hopping operators which we need to 

keep track of, ^ and / * , . 

The initial block states and the matrix elements of the hopping operators can be constructed 
immediately. A Mathematica file, kondo_Ul_charge_Ul_spin.nb, with all the numerical details 
of this calculation as well as an input file is provided by the package. The corresponding block 
and the local states are listed in Table l/TTl 

In the first iteration, the reduced matrix elements of the Hamiltonian matrix between the block 
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states (Tabic I2.1[) are given by 
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/ n 


n 
u 


n 
u 


n 
u 


n 

u 


u 


u 


n \ 
1 


U 


u 


U 


U 


U 


n 
u 


U 


U 


o 


o 


1/4 


() 


() 


() 


() 


o 











-3/4 


























1/4 


























1/4 
































V o 




















o J 



(2.10) 



In general there might be off-diagonal matrix elements in this matrix. The code takes care of 
this and diagonalizes o(u||iJo||w)o before proceeding to the next iteration. 
The matrix elements of the /q ♦ are 



oHI/o,t||«}o = 



while those of /q , are given as 
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(2.11) 



(2.12) 



These matrix elements are needed to construct the hopping terms in the first iteration. In 



addition, one also needs the matrix elements of the hopping operators j 
the local states, listed in Table 12.11 These are given by 
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(2.13) 



Finally, if we want to determine the auto correlation function of S z , e.g., then we also need to 
compute the matrix elements this operator, 
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(2.14) 
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Table 2.2: Block states (left) and local/added states (right) for the single channel Kondo 
model when J7 c ^ arffe (l) x SU sp i„(2) symmetry is used. The first block states are formed from 
the impurity spin and the conduction electron at site n = of the Wilson chain. 



We remark at this point that, while S z is an irreducible tensor operator and it has a spin 
quantum number 0, the operators S x and S y are not irreducible tensor operators with respect 
to the UgpiniX) spin rotations, since they transform among each other. They are, however, 
linear combinations of the irreducible tensor operators, S±, which have spin quantum numbers 
±1 under rotations around the z-axis. 

Once all this information typed into the input file, input.dat, the flexible DM-NRG code is 
ready to compute the finite size spectrum, the spin's spectral function, and the real part of 
the retarded spin-spin correlation function using both traditional NRG as well as DM-NRG 
methods! An example input file for this case is provided in the directory input_f iles. 



Kondo model with U c harge{l) x SU sp i n (2) symmetry 



Let us now show on the example of the Kondo model, how non-Abelian symmetries can be 
used. Let us thus enumerate all necessary information that is needed to perform a calculation 
for the spin spectral function of the single channel Kondo model with V c harge[l) x SU sp i n {2) 
symmetry. A Mathematica file (kondo_Ul_ch.arge_SU2_spin.nb) where these inputs have been 
computed is also provided by the package. 

In the first iteration we have four different block multiplets, u = 1,...,4 formed from the 
impurity spin and the conduction electron at site n = 0, while there are three added multiplets 
in each iteration, \i = 1,2,3. These states and their quantum numbers are listed in Table 
12.21 (only highest weight states are given). The reduced matrix elements of the Hamiltonian 
between the block states then read 
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(2.15) 



In this case, the matrix elements of /q ^ and /q , are related with each other by symmetry, and 
therefore they form a single hopping operator multiplet /q of spin 1/2. The matrix elements of 
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Table 2.3: Block states (left) and local/added states (right) for the single channel Kondo model 
when SU c harge(2) x SU sp i n (2) symmetry is used. The first block states are formed from the 
impurity spin and the conduction electron at site n = of the Wilson chain. 



/g need to be determined using the Wigner-Eckart theorem, and are given by 
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(2.16) 



To generate the hopping, we also need to know the reduced matrix elements of /],[ between 
the added (local) states. These are computed as 




(»\\fl+i\W) ( 1 

-y/2 



(2.17) 



Finally, to compute the spin-spin correlation function, we need to compute the reduced matrix 
elements of the impurity spin. As discussed before, the components of the spin operator are 
grouped in this case into a single tensor operator, S. The reduced matrix elements of this 
operator are computed as 



o(u||S||u)o = 



/3/i 

-V3?4 
1/2 l/y/2 

ooo vW4 / 



(2.18) 



To run the code with U c harge(X) x SU sp i n (2) symmetry, one only needs to type this information 
into the input file, input.dat, and then run the code. You can find the corresponding input 
file, kondo_model_Ul_c_SU2_s.dat in the directory input Jiles. 
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Kondo model with SU charged) x SU sp i n (2) symmetry 

The Hamiltonian (|2.2[) is invariant under the action of SU c harge(2) rotations in charge space, 
Uc = e iZ Q Q : generated by the operators Q x = (Q+ + Q~)/2, Q y = (Q+ - Q~)/2i and Q z with 



»=o V={T,i} / 



(2.19) 



(Q 



and parametrized by real, three-component vectors luq. Since the spin symmetry generators 
commute with the charge symmetry generators, the Hamiltonian (|2.2p possesses a symmetry, 

SU c harge{2) X SU S pi n (2) . 

Charge SU(2) symmetries must always be used with care, and one must always carefully check if 
a given set of operators forms an irreducible tensor operator. One can check, e.g., that the cor- 
rectly defined hopping operators read, 7* = {{ln)ar} = {fl A , fl d , (-1)™ f n ,l> (~1)™ +1 /n,t}> 
and that they indeed transform as a set of spin 5=1/2 and charge spin Q = 1/2 irreducible 
tensor operators. Notice the curious sign-dependence of the last two operators and their relative 
signs. Similar factors of (—1)™ appear in the states generated from the highest weight states in 
Table l2~3l through the action of the operator Q~ . 

In the first iteration we have three different block multiplets, u = 1, ... ,3 formed from the 
impurity spin and the conduction electron at site n = 0. There are two added (local) multiplets 
in each iteration, fx = 1,2. These states and their quantum numbers are listed in Table [2~31 (only 
highest weight states are given) . The reduced matrix elements of the Hamiltonian between the 
block states then read 

/ -3/4 \ 
o(u\\H \\v) = J [ 1/4 . (2.20) 

\ 00/ 

The matrix elements of need to be determined using the Wiegner-Eckart theorem, and are 
given by 

/ -V2 



V V2 V 2 



-V2 




(2.21) 



To generate the hopping, we also need to know the reduced matrix elements of 7} between the 
added (local) states. These arc computed as 

- ( ^ f ) ■ (2-22) 

Finally, to determine the spin spectral function, we need to compute the reduced matrix ele- 
ments of the impurity spin operator: 

/ -VW^ \ 
oH|5||«) = 1/2 1/V2 . (2.23) 
V / 

To use the code with SU c harge(2) x SU sp i n (2) symmetry, one only needs to type this information 
into the input file, input . dat, and then run the code. You can find the corresponding input file, 
kondo_model_SU2_c_SU2_s.dat in the directory input_f iles, and a Mathematica notebook to 
compute these matrix elements is also provided. 
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2.1.2 Extension to arbitrary symmetries 

The procedure outlined above can be generalized to essentially any number and type of non- 
Abelian symmetries and to any Hamiltonian of the form: 

oo 

+ «n+i) • (2-24) 

n=0 

The first term in Eq. (|2.24|) . H.q, describes the quantum impurity coupled to the 'environment'. 
This is thus the 'interaction part' of the Hamiltonian. The dynamics of the environment is 
described by hopping terms that connect nearest-neighbors only, T n<n +i, and the on-site terms 
Ti. n . As already stated before, these on-site terms can account for the breaking of electron-hole 
symmetry or for the presence of superconducting correlations [T5] , 

The above Hamiltonian can also be diagonalized iteratively using the recursion relation, (|2.4[) . 
The procedure is essentially the same as the one outlined in the previous subsection, with the 
modification that we now may have V different symmetries, which commute with the Hamil- 
tonian and each other, and which we shall label by an index 7 = 1, . . . ,T. Correspondingly, 
multiplets will be labeled by T different quantum numbers, {Q 7 } = {Q 1 , ■ ■ ■ , Q r }, with Q 1 the 
label of the irreducible representation of symmetry 7. States within a degenerate multiplet can 
then be labeled by internal quantum numbers, which, by analogy to the spin operators, we shall 
denote by QJ. Finally, the representation indices and the internal quantum numbers can both 
be grouped into vectors, Q_ = {Q 1 , Q 2 , Q r } and Q_ z = {Q\, Q 2 Z , Q T Z }. 
Similar to multiplets, a set of irreducible tensor operators, {A} is also characterized by some 
quantum numbers a, and a z . The hopping terms, r nj „ + i, can also be decomposed using a set 
of such irreducible operators, which we shall call 'hopping operators', /' — > C, and label them 
by A. Thus the hopping operator multiplet, C\ has quantum numbers c A and members of this 
multiplet are labeled by c\. Ordering the members of a hopping operator multiplet into vectors, 
= {C CXiC *}, the hopping can be simply written as 

r„,„+i = ' £l+i,A + h A ■ (2-25) 



Note that here different hopping operator multiplets may have different hopping amplitudes, 
and also, remember that now the C_ n A 's denote creation operators. An example for having two 
operators has already been for the Kondo model with t/ c /j ar9e (l) x U sp i n (l) symmetries, but in 
most cases several hopping operators must be used for multichannel models too. 
In our general framework the iteration step n — * n + 1 reads now as follows: 

• Previous block states u can be characterized by their symmetry indices and labels, Q , Q z , 
u,Q u ,Q z u ) . Similarly, local states ji are also characterized by their symmetry indices 



% = and intcrnal labels gl = {<il,z}- 



The new basis states can be constructed using Clebsch-Gordan coefficients as 



where we have also introduced a general notation for the Clebsch-Gordan coefficients as 
a product of coefficients, one for each symmetry in the problem 



q q z ;Q Q 2 

—fA—fJ, ——It— —I 



r 

Q&) = II {il^QlQl, z \QlQh) ■ ( 2 - 27 ) 

7 = 1 
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• Knowing the reduced matrix elements, n (u |j C n .\ \\ v) n and (/i || C n +\,\ \\ v), and the 
eigenvalues E™ from iteration n one can construct the Hamiltonian in the new basis. 

• Diagonalizing this Hamiltonian, one can determine the new eigenstates, u, Q^,Qt 

— u — u/ n+1 

their energies, E™ +1 , and the irreducible matrix elements, n +i(u || C„+i,a || v)n+i- If a 
spectral function is computed then the corresponding operator's reduced matrix element, 
n +i(u || A || must also be computed. 

• After the iteration some of the high-energy states are also discarded before one proceeds 
to the next iteration. 

Our flexible DM-NRG code implements this rather general procedure dynamically, and in such 
a way that various symmetries can be taught to the code. If you want to teach a new symmetry 
to the code, feel free to contact us for help. 



2.2 The DM-NRG method 

Wilson's procedure gives accurate results for the finite size spectrum, and it can also be used 
to compute spectral functions [13] . However, there are several problems with the simple ex- 
tension of Wilson's method. On one hand, in the original approach, spectral functions have 
been computed using only the approximate ground states in a given iteration. This may give 
incorrect results if the fine structure of the ground state shows up at subsequent iterations only. 
Hofstetter proposed to use a density matrix NRG (DM-NRG) method to treat this problem 
[TJ. Furthermore, in the original method one simply chops off parts of the Hilbert space, and 
therefore spectral sum rules are not respected. This problem has been solved only recently, 
through the introduction of a complete basis set [5]. In this approach, rather than constructing 
approximate eigenstates for a chain of length n, one constructs approximate eigenstates of H n 
that live on a chain of length TV > n, by adding environment states 

u 7 Q ,Q* 



u,Q Q z ;e) (2.28) 
where e labels all possible states of the environment, i.e., chain sites from site n + 1 to N. 
The reduced density matrix 

Physical quantities are then computed from the density matrix, which is approximated as 

JV 



Q = (2-29) 
^E E E^" ^9.>)_UQ a ,Ql;e , (2.30) 



z 

e u^Discardcd Q z 

where the partition function is defined as 

Z = Y^ dim£-" dim(u) e^" , (2.31) 

u 

where dim(u) is the dimension of the multiplet it, and dimi oc = dim(/z) denotes the dimen- 
sion of the local basis at a given site of the chain. The second summation in Eq. (|2.31|) runs 
over states that have been discarded in iteration n — > n + 1 . 
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The calculation of a spectral function then consists of a forward sweep, where we determine all 
cigenstatcs in a given iteration, and a backward sweep, where we determine g and the truncated 
reduced density matrices, defined as 

flM = Tv N _ n { Qm) ■ (2.32) 
The details of the corresponding recursion relations were given in Ref . [5] . 



Spectral function calculations in the fiexible-DM-NRG framework 

The main purpose of the code is to compute the retarded Green's function of two local operators, 
A and B. This is defined for two irreducible tensor operators as 



G 



(*) = -i ( [Aa,a z (t), 4^ (0)] } 6 (t) . (2.33) 



where the quantum numbers a, a z and 6, b z refer to the quantum numbers of operators and 
. For fermionic operators one must replace the commutator above by an anticommutator. By 
symmetry, this Green's function is non-zero only if the spectral operators A and B transform 
accordingly to the same representation, 



G n A flt (t) = G2,Bt(t)Sa,b^, K ■ (2.34) 

Although it is a cumbersome expression, for completeness, let us give here the expression 
obtained generalizing the procedure of Ref. [5] ■ 

g% b m = Y. E E «<<n* w iii> B 

n=0 ieD,K (j.fe)^(K,K) 



z+ \{E2 - El dim(o) z — \{Ef + E^) — E% dim(o) 

(2.35) 

Remarkably, this formula contains exclusively the reduced matrix elements and the dimensions 
of the various multiplcts and operator multiplets. Here the second sum is over all the multiplcts 
k of the given iteration subject to the restriction that j, k do not belong to kept states at the 
same time. No summation is needed for states within the multiplcts. The sign factor £ is £ = 1 
for bosonic operators, while it is £ = — 1 for fermionic operators, and dim(a) = Il^i dim(a 7 ) 
is the dimension of the operator multiplet A\ a . The code computes the spectral function 
using this formula, and the full Green's function is then computed by doing a Hilbert transform 
numerically. 

Smoothing procedure 

In general, the spectral functions are given as weighted sums of 5-functions of the form: 

A AB t(uj) = --3mG^ t (w + «0+) (2.36) 
= 5{u-u)i) (2.37) 

i 

where the weighting coefficients Wi can be expressed from Eq. (|2.35| after analytical continu- 
ation to the real axis. To get a smooth spectral function, however, the delta functions need to 
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be replaced by some broadening kernels K{u>,u>i) that fulfill the condition 



du K(u, uJi) = 1 . (2.38) 

This condition guarantees that spectral sum rules remain valid after the broadening procedure. 
In our code we use two types of kernels to generate the real axis spectral function: (a) the 
log-Gaussian kernel [16] and (b) the interpolative kernel that we constructed to avoid certain 
difficulties with the interpolation around lo = 0. 
For the log-Gaussian method the following kernel is used 

K log (u>, Ui ) = 1 e -» a /4 e -0na,-ln«*)V&»l. ( 2 . 3 9) 

This kernel works rather well at T = calculations, however, it extrapolates to the u> — > ±0 
values in a singular way. Therefore, it has problems at finite temperatures and also in cases 
where spectral losses lead to an artificial jump in the spectral function at lo = 0. 
The interpolative scheme avoids this difficulty by using a kernel 

K mt (^) = J- e -[*M-(-)]V^ (2 .4 0) 

with the smoothening function defined as 

1 / \ / / \ 2 



,H = -tanh^-j ln^-j + e-J , (2.41) 

by choosing 7 w Tq e.g. This formula interpolates smoothly along the real axis for all fre- 
quencies. It has the property that for frequencies smaller than the 'quantum temperature', 
\lu\ <C Tq, the broadening is Gaussian, while for \lo\ 3> Tq it becomes log-Gaussian. In this 
scheme, for zero temperature calculations the quantum temperature Tq must be the smallest 
energy scale in the problem, on the other hand for finite temperature calculations, the quantum 
temperature must be set to be in the range of the temperature itself, Tq ~ 0.5T -j- T. In both 
interpolative methods there is a fitting parameter b, which controls the widths of kernels and in 
general is A dependent, which should be in the range of b ~ \/A/2 (for A = 2, b = 0.5 -j- 0.9). 
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Chapter 3 

Installation and technical support 



3.1 Memory and disk space requirements 

The resources required by the DM-NRG code depend a lot on the type of the ongoing calculation. 
Increasing the number of block states that are kept during a run, the CPU time needed for 
performing the calculation, the necessary memory, as well as the used disk space may increase 
substantially. As an example in Table [5TTI we present some estimates for a simple model. The 
following information was extracted by running the code on an Intel(R) Core(TM)2 Duo CPU 
T7250 @ 2.00GHz processor with 4MB CPU cache. We have used the CentOS 5.1 distribution 
for the Linux operating system. The code was run by using an input file for the Anderson 
model with U(l) x SU(2) symmetries that comes with the distribution. We made a run for a 
flat band and performed 50 iterations. The compiler used was gec version 4.1.2. 



Model name 


number of 


CPU-type 


CPU time 


disk space 


memory 




block states 




(sec) 


(MB) 


(MB) 


Anderson model 


50 


Intel Corc(TM)2 Duo 


1.27 


12 


4 


U(l) x SU(2) 


100 


CPU T7250 @ 2.00GHz 


3.56 


25 


7 


symmetries 


200 


with 4MB cache 


14.39 


69 


12 




500 




156.26 


341 
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Table 3.1: CPU time, memory and disk space requirements. 



Important note: when temporary files are also generated in ASCII format by setting text_swap_files_flag 
= ON in <SECTION-FLAGS> of the input file input .dat, approximately four times larger 
disk space is required than the ones shown in the table. 

3.2 Compilation environment 

The code can be compiled with either g++ or with ice (intel C++ compiler). The Intel com- 
piler, ice is available for free for non-commercial and non-academic use after signing a license 
agreement. 

3.3 Libraries 

The Flexible-DM-NRG code is linked with state of the art libraries. The following list of 
libraries (packages) are needed to have the f lexible-dmnrg package properly installed. The 
easiest way to install these packages is by using the rpm - RPM package manager method. 
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For further information and links to the websites from where these libraries can be downloaded 
please refer to the web page of the code http: //www. phy.bme.hu/^dmnrg/index. html. 
List of required libraries: 

• libg2c - Fortran 77 Libraries 

• bias - Basic Linear Algebra Subprograms 

• lapack - The LAPACK libraries for numerical linear algebra. 

• gmp - A GNU arbitrary precision library 

• gsl - GNU Scientific Library 

• gsl-devel - GNU Scientific Library - development files 

3.4 Automatic installation 

Once the libraries were successfully installed, run the following shell script in the main directory 
of the code: 
. / configure 

This command must run without errors. It checks your system, searches for the necessary 
libraries, and sets the library paths. An error message returned while running ./configure 
means that some packages might be missing or not properly installed. If the ./configure has 
run without errors then the make . sys and make . rules and conf ig.h files have been generated. 
These files now contain the system variables needed for the installation 0. 
Next run the following command: 
make all 

If everything goes smooth and no error messages were generated, then the code is successfully 
installed. If for some reasons the automatic configuration does not work, please, follow the 
instructions in the next subsection and install the code manually. If you do not succeed, feel 
free to contact us. 

3.5 Manual installation 

Sometimes, for some mysterious reasons, the autoconf does not work properly and some of the 
libraries are not detected. In this case, the make . sys file, that can be found in the main folder 
of the code, needs to be edited manually. Altogether there are ten shell variables that need to 
be edited. 

• CXX - the name of the C++ compiler, is usually g++ or ice. 

• C - the name of the C compiler, usually gcc. 

• CC_FLAGS - flags for the g++ or ice compiler (this variable can be left empty at this stage). 
You can set here the optimization level to speed up the code. 

• GSL_FLAGS - path to the gsl include files, usually -I/usr/include. 

• GMP_FLAGS - path to the gmp include files, usually -I/usr/include. 

• FLIBS - fortran libraries, usually g2c or gf ortran. 

• LAPACK_LIBS - lapack library, usually -llapack. For static libraries the full path name needs 
to be provided sometimes: /lib/liblapack. a. 

• BLAS_LIBS - bias library, usually -lblas. For static libraries the full path name needs to be 
provided sometimes: /lib/libblas . a. 

• GSL_LIBS - gsl library, usually -lgsl -lgslcblas. For static libraries the full path name 
needs to be provided sometimes: /lib/libgsl . a. 

• GMP_LIBS - gmp library, usually -lgmp. For static libraries the full path name needs to be 
provided sometimes: /lib/libgmp . a. 

lr The automatic installation checks for the g++ and gcc compilers and does not check for the ice. To run the 
code with ice, manual intervention is necessary (see section Manual Installation), and the CXX environmental 
variable needs to be modified from g++ to ice. 
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Once the make . sys file is properly edited run 
make all 

If the compilation runs without errors then the code is successfully installed and all the binaries 
have been generated. 

3.6 Generated binaries 

After a successful installation the following binaries are generated: 

• fnrg - the main program that does the NRG/DM-NRG calculations. 

• sfb - utility that does the spectral function broadening. 

• he - utility that computes the hopping for an energy-dependent density of states. 

• es - utility to generate/analyze the energy spectrum. 

• cgc - utility to generate/compute the Clcbsch-Gordan coefficients 

To run the code the input.dat file must be set properly. For further details, refer to the 
next sections of the manual. Example input files are provided with the code for the single 
band Anderson model with U(l) charge and spin symmetries and for the Kondo model with 
combinations of U(l) and SU{2) symmetries. These 'ready to run' input files are available in the 
. /input Jriles folder. Furthermore, we also provided the Mathematica files which were used 
to compute the content of these input files. Some of the technical details of these calculations 
were summarized in Chapter [21 

3.7 BUG reports 

We have made this code publicly available, so if you think you have found a bug in the code, 
please let us know so we can fix it for future releases. 

Please send any bug report to: dmnrg@neumann.phy.bme.hu. Use BUG REPORT in the 

subject area. 

Please include the followings in any report: 

• The version number of the code. 

• input.dat file so we can test the code with it. 

• A description of what is wrong. If the results are incorrect, in what way. If you get a crash, 
etc. 

• Please do not send core dumps, of executables. 

• The name of the compiler that was used and its version. 

• The output from running ./configure 

• If the bug is related to configuration, then the contents of config.log. 

If the bug report is good than we will do our best to fix it. If the bug report is poor we will 
not do much about it, probably just ask for a better report. 

If you think something in this manual is unclear, or downright incorrect, or if the language 
needs to be improved, please send us a note to the same address. 

3.8 Some hints for developers 
3.8.1 Organization of the code 

The components of the program are located in three main directories. Directory . / src/ contains 
the *.cc source files, ./include/ the corresponding .h files and . /input .files/ the user 



23 



defined input files. The directory ./doc/ contains the documentation of the code in LATEX 
and HTML forms. Table I3T21 gives a hint on the structure of the code. 



Makefile 
DDXYFILE 

sre/ *.cc 
include/ *.h 
input Jiles/ *.dat 
results/ 

doc/ nrg_user_guide.tex 
html/ 
tex/ 



configuration for on-line documcntion 
source code (source files) 
source code (header files) 
directory for input files 
output directory of the code 
documentation of the code 
on-line documention 
on-line documcntion 



Table 3.2: Directory structure of the code. 

The other files in the main directory are configurations files needed for properly installing the 
code. 

3.8.2 On-line documentation of the code 

The documentation of the source code is managed with the use of the doxygen a free software 
under Linux. At present version 1.4.6 is used. The configuration of doxygene is set in file 
Doxyfile and the program is called from shell as doxywizard. Description of the classes, variables 
and functions is included directly in the source code (in the *.cc and *.h files), therefore, all 
comments about new variables, functions, classes etc. are documented right away. The on-line 
documentation file ./doc/html/. . . .html can be opened with any WEB-browser. The DM- 
NRG reference manual in pdf form (refman.pdf) can be generated by the command: latex 
. / doc/latex/ ref man . tex. 

All the information about the functions (classes) i.e., what it does, how to call, description of 
the input and output variables must be located in the main body of the functions, otherwise, 
this information will not be accessible in the output files generated by doxygene! The hierarchy 
of the classes, dependencies and calling functions are generated automatically by doxygen. 

3.8.3 Version management 

Version management of the code is carried out with the use of Subversion system(svn) a 
free software under Linux. This program is similar to Control Version System (cvs) but 
has several additional features. The WEB page of the svn interface can be accessed by the 
members of the flexible DM- NRG developers group only. 

The management of this manual - in LATEX form-, however, has to be done manually. There are 
no scripts yet, that would automatically update the manual once the code has been changed. 
However, since all directories and files - including this manual as well - arc checked in to the 
svn it is straightforward to update the manual and all changes of it are documented by the svn 
itself. 
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Chapter 4 

Using the DM- NRG code 



4.1 Main features of the code 

If the installation was successful then, to run the code for your problem of interest, you first 
need an input file, input.dat. The structure and content of this file shall be described in 
detail later in this chapter, and a few examples are also provided with the code in the library, 
. /input Jiles. Furthermore, we provided some Mathematica files in the mathematical iles 
folder, which were used to generate these input files. For the generation of the irreducible 
matrix elements it is sometimes important to know, what is the definition of the Clebsch- 
Gordan coefficients used by the code. For this purpose, a utility called cgc is provided in the 
package that computes the Clebsch-Gordan coefficients. 

When running the code, all the results are written into the folder ./results. The content of 
these files shall also be described later in this chapter. Finally, these results can further be 
analyzed with the utilities sfb - real axis spectral function generator and es - energy spectrum 
generator, also provided with the code, and shall be described later. 

Main features included in version 1.0.0: 

• Zero and finite temperature calculations of the spectral function of not necessarily the 
same operators. 

• Static averages. 

• NRG type calculations or DM-NRG computations (using the complete basis set). 

• The code uses a single input file that can be easily modified. Examples are available for 
Anderson as well as Kondo Hamiltonians in the presence of U(l) and SU(2) symmetries 
for spin as well as charge symmetries. 

• Flat density of states as well as arbitrary energy-dependent density of states can be con- 
sidered. In the latter case arbitrary precision routines are used to generate the parameters 
of the Wilson chain. 

• Utility for broadening the spectral function is provided to compute the imaginary and the 
real parts of the retarded Green's function 

G% B ,(t) = -i([A(t), flt(0)] ± )e(i), (4-1) 
with the operators A and B not necessarily identical. 
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• Utility for computing the on-site energies as well as the hoppings for a non-uniform density 
of states is provided. 

• Utility to generate/analyze the finite size spectrum is provided. 

Utilities provided in the package: 

• he - generates the hopping amplitudes along the chain and the on-site energies when an 
energy-dependent density of states is used. (This utility provides input files that are read 
directly in the fnrg binary code). 

• sfb - does the broadening of the spectral function and computes the real part of the 
Green's function by a Hilbert transform. This utility uses the outputs of the fnrg to 
generate these data. 

• notebook examples - there are some extra Mathematica 6.0 notebook examples where 
the matrices for Hamiltonians as well as for some hopping and spectral operators are 
computed. These notebooks can easily be extended to other impurity models. 

• es - generates the energy spectrum. This utility uses the output of the fnrg to generate 
the data. 

• cgc - an interactive utility, which generates the Clebsch-Gordan coefficients used by the 
code. 



4.2 Running the code 

In this section we shall give some information on how exactly to run the code. At this point we 
suppose that the code was successfully compiled and linked properly with the required libraries. 
If this is the case there should be five binary codes in the main folder, fnrg, sfb, he, cgc, and 
es. Depending on the type of calculation that is performed, the following steps need to be 
followed: 

4.2.1 Flat density of states 

Most NRG calculations are done for a flat density of states, g(co) = |, with ui G [—1,1]. In this 
case there is no need for the he routine, since the hoppings along the chain are computed on 
the fly. In this case the on-site energies are all zero as a consequence of electron-hole symmetry. 
For this type of calculation, first the hoppings-onsite-energies-flag in the <SECTION-FLAGS> 
needs to be set to 'OFF' in the file input.dat. Then the fnrg must be run without any 
arguments in the command line. The name of the input file must always be input.dat, since 
the code looks for this particular file to read in the data. 

Outputs are saved in a folder in results/. Results for the spectral function can be generated 
by running the sfb binary without any argument after exiting normally from the run. Note 
that sfb reads the same input file, input.dat, so do not modify it between the two runs 
except for the part <SECTION-SPECTRAL_FUNCTION_BROADENING>. Depending on 
the flags set in <SECTION-SPECTRAL_FUNCTION_BROADENING>, spectral functions, 
Green's functions and static averages are computed. Note that the real part of a Green's 
function can only be computed if the corresponding spectral function has previously been 
computed. 
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4.2.2 Energy dependent density of states 

The user can also provide an arbitrary density of states (DOS) for the code. In the present 
version of the code, the DOS is read from a file, dos.dat. This file must be in the folder, 
. /dos_mapping, and it must be provided by the user. This file must contain at least two 
columns, one of them being the energy, and the others being the corresponding density of 
states. There should be as many density of states columns as hopping operators are defined in 
the input file. 

There are some constraints that the DOS must satisfy 

• The DOS needs to be rescaled such that all energies are in the interval, J? 6 [-1,1]. In 
other words, the energy unit of the NRG calculation (=1) must be larger than the range 
of the density of states. 

• The DOS needs to be normalized to 1. 

• The energies must be ordered. 

• Each density of states column corresponds to a given hopping operator in the input . dat 
file. The columns must be ordered accordingly to the order of the hopping operators in 
the input.dat file. 

There is no further requirement for the mesh, both uniform and non-uniform meshes can be 
provided. 

When the NRG calculation is performed with an energy-dependent DOS, the hoppings and 

the on-site energies must be computed first. This must be done interactively, using the binary 

he. When running he, one needs to specify an iteration number. This can be larger than 

the number of iterations specified in the input . dat file. After a successful run of he, the 

hopping amplitudes and the on-site energies are saved in files having descriptive names such as 

hopping_couplings_* . dat and on_site-energies-*.dat in the folder . /dos_mapping/. 

Next, the input . dat file needs to be set-up: There, make sure that the hoppings_on_site_energies jf lag 

in the <SECTION-FLAGS> is set to 'ON' for this type of run. 

Once the hopping amplitudes and the on-site energies were computed, the same steps must be 
followed as in the case of a flat density of states. The code fnrg must be run without arguments. 

Note : The routine he is using arbitrary precision libraries for computing the hoppings and 
the on-site energies, however due to the accumulation of the errors during the calculations 
sometimes the he utility may give unexpected results for the hoppings. Therefore, we advise 
the user to check the results for the hoppings before pluging them into fnrg for any kind of nrg 
calculations. 



4.3 Initialization and the input file 

In the initialization step - the zero'th NRG iteration step - one has to define the Hamiltonian. 
Therefore, for the iterative process, among others, the user has to provide the following input 
parameters: 

1. One needs to specify the number (symmetry _no) and type of symmetries. 

2. One needs to specify the number (block_statejnumber) of initial block states, i.e., multiplets 
that span the Hilbert space for n = (we denote this number in this manual by K) . One 

also needs to define the quantum numbers Q of the block states, u, Q ,Q Z ) ■ 

U U U j q 

3. Block Hamiltonian matrix - The matrix elements of the K x K Hamiltonian matrix Hq, 
defined on the K block states. 
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4. block hopping operators - The reduced matrix elements of all hopping operators between 
the block states. These are K x K matrices. 

5. spectral operators - The reduced matrix elements of the operators, whose spectral functions 
we compute. These are also K x K matrices. 

6. local state number - The number of multiplets that span the local Hilbcrt space (we label 
this number L). 



8. local hopping operators - The reduced matrix elements of the hopping operators. These 
arc L x L matrices. 

This list is rather incomplete, and explains only briefly the meaning of the most basic vari- 
ables. In addition, many other, less relevant parameters must be specified for the NRG process 
(broadening, maximum iteration number, etc.). All these input parameters are read from file, 
input.dat, which must be constructed for every model separately. Some examples of such 
input files can be found in the folder, . /input_f iles, and also in the appendix of this manual. 
A detailed description of these parameters shall be provided in forthcoming subsections. 

Structure of the input file 

The input file is organized along sections, and contains a lot of comments. As a general rule, 
any line that starts with the symbol is treated as a comment line in these input files and 
is ignored by the code. You can thus add additional comment lines for yourself using this 
convention. 

The input file has the following overall structure (that we shall describe in detail in subsequent 
subsections) : 



1. 


<SECTION- 


PARAMETERS 


2. 


<SECTION- 


FLAGS > 


3. 


<SECTION- 


SYMMETRIES > 


4. 


<SECTION- 


BLOCK_STATES> 


5. 


<SECTION- 


LOCAL_STATES> 


G. 


<SECTION- 


LOCAL_STATES.SIGNS> 


7. 


< SECTION- 


BLOCK_HAMILTONIAN> 



(a) <BLOCK_HAMILTONIAN_TERM> 

i. <BLOCK_HAMILTONIAN_NAME> 

ii. <BLOCK_HAMILTONIAN_COUPLING> 

iii. <BLOCK_HAMILTONIAN_REPRESENTATIONJNDEX> 

iv. <BLOCK_HAMILTONIAN_MATRIX> 

8. <SECTION-HOPPING.OPERATORS> 

(a) <HOPPING.OPERATOR> 

i. <HOPPING_OPERATOR_NAME> 

ii. <HOPPING.OPERATOR_REPRESENTATIONJNDEX> 

iii. <HOPPING.OPERATOR_SIGN> 
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local states - Quantum numbers q of the local states, /i, q ,q z 
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iv. <HOPPING.OPERATOR_MATRIX> 



9. <SECTION-SPECTRAL_OPERATORS> 

(a) <SPECTRAL_OPERATOR> 

i. <SPECTRAL_OPERATOR_NAME> 

ii. <SPECTRAL_OPERATOR_REPRESENTATIONJNDEX> 

iii. <SPECTRAL_OPERATOR.SIGN> 

iv. <SPECTRAL_OPERATORJVIATRLX> 

10. <SECTION-STATIC_OPERATORS> 

(a) <STATIC.OPERATOR> 

i. <STATIC_OPERATOR_NAME> 

ii. <STATIC_OPERATOR_REPRESENTATION_INDEX> 

iii. <STATIC.OPERATOR_SIGN> 

iv. <STATIC.OPERATORJVIATRIX> 

11. <SECTION-SPECTRAL_FUNCTION> 

12. <SECTION-SPECTRAL_FUNCTION_BROADENING> 

13. <SECTION-LOCAL_HOPPING_OPERATORS> 

(a) <LOCAL_HOPPING.OPERATOR> 

i. <LOCAL_HOPPING_OPERATOR_NAME> 

ii. <LOCAL_HOPPING_OPERATOR_REPRESENTATIONJNDEX> 

iii. <LOCAL_HOPPING.OPERATOR.SIGN> 

iv. <LOCAL_HOPPING.OPERATOR.MATRIX> 

v. <LOCAL_ON_SITE_ENERGY_MATRIX> 

4.3.1 <SECTION-PARAMETERS> 

In this section, the main parameters of the given model and also those of the DM-NRG method 
are set. The order of the parameters does not matter as long as none of them is missing. 

Model parameters 

• model: The name of the model under investigation. It can be any string. This name will be 
automatically included in the name of the folder, where output data are saved. 

• lambda: The value for the logarithmic discretization parameter, A. It is usually fixed between 
2 and 3. 

• symmetry .no: The number of the symmetries used. The present version of the code handles 
U(l), 5/7(2), Z(2), chargeSU(2) symmetries. The properties of the symmetries are set in 
< SECTION-SYMMETRIES > by defining their names and setting upper and lower bounds on 
their representation indices. Later versions of the code will include SU(3) symmetry and crystal 
field symmetries too. 

• coupling_no: The number of independent coupling constants. In general, the local Hamil- 
tonian Hq can contain several couplings. The on-site part of the Anderson model, e.g., can be 
written as 

H = U(n - l) 2 + (e d + U/2) n + V + h.c.) , 
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i.e. it contains three different couplings, U , Cd and V. These couple to the operators, (n — l) 2 , n, 
and (tfjr/o.CT + ^-c). Thus the local Hamiltonian, Hq, can be written as a sum of coupling-no = 3 
independent terms, 

coupling-no 

H = J2 Gs^s, (4.2) 

5=1 

where Tts's denote the independent local terms and Gg stand for the corresponding couplings. 
The matrix elements of H s are specified in <SECTION-BLOCK_HAMILTONIAN>. 
Important note: The number of <BLOCK_HAMILTONIAN_TERM> items must agree with 
the 'coupling_no' set in this section. 

• spectral_operator_no: The number of the operators for which one wants to compute 
the spectral functions. The detailed parametrization of these operators is carried out in 
<SECTION-SPECTRAL_OPERATORS>. The number of <SPECTRAL_OPERATOR> items 
must agree with the 'spectral_operator_no' set in this section. 

• static_operator_no: The number of operators for which static averages are computed. The 
detailed parametrization of these operators is carried out in <SECTION-STATIC_OPERATORS>. 
The number of <STATIC_OPERATOR> items there must agree with the 'static_operator_no' 
set in this section. 

• hopping_operator_no: The number of independent hopping operators. The detailed parametriza- 
tion of these operators is carried out in <SECTION-HOPPING_OPERATORS> and <SECTION- 
LOCAL_HOPPING_OPERATORS> . 

• block_state_no: The number of multiplcts that span the Hilbert space for iteration n = 0. 
The detailed parametrization of the states is carried out in <SECTION-BLOCK_STATES> 
by setting the representation indices of all the states corresponding to all the symmetries used 
during the calculation. 

• local_state_no: The number of local multiplcts, added in each iteration. In the present 
version of the code this is the same for all sites, < i < N. The detailed parametrization of 
local states is carried out in <SECTION-LOCAL_STATES>. 

• local_coupling_no: The number of independent terms in the on-site Hamiltonian operators, 
Tt„. In the presence of electron-hole symmetry of the conduction band and without any other 
correlations (superconducting like) the on-site terms are all zero. If we have only on-site energies, 
£„, then local_coupling_no= 0. 

• spectral_function_no: The number of the spectral functions that should be calculated. 
One can compute several spectral functions within the same run, at the expense of using 
more memory. The detailed parametrization of the operators is carried out in <SECTION- 
SPECTRAL.OPERATORS> . 

DM-NRG parameters 

• max_state_no: The upper limit for the number of kept block multiplets. 

• iteration_no: The maximum number of the NRG iteration steps, N . The program terminates 
once this value is reached. 

• intervaLno: The energy grid used in each NRG iteration step when spectral functions are 
calculated. Note that in DM-NRG a very large number of delta peaks is generated in every 
iteration. Therefore, these delta peaks are stored on a finite mesh. In the present version of 
the code, the difference between the maximum and minimum of the energy spectrum at each 
iteration is divided by interval-no to obtain the size of the mesh, du>. Typically interval-.no ~ 
1000. 

• degeneracy_threshhold: If two states are such that the difference between their energies 
is less than degeneracyjthrestLhold these states are considered to be degenerate. 

When discarding states, those lowest in energy are kept while the higher energy states are 
discarded, such that at the next iteration no more than max_state_no are kept. Degenerate 
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states are block discarded. 

• temperature: The temperature used for the calculation. When the temperature is fixed 
to zero, temperature = 0, the code runs by using the number of iterations provided by the 
user. In case of a finite temperature, temperature > 0, the variable iteration-no is determined 
and modified during the run to correspond to the temperature specified, using the well-known 
relationship, T n ~ (1 + ±) A~"/ 2 . 

4.3.2 <SECTION-FLAGS> 

In this section various flags controlling the DM-NRG method are set. 

• dmnrg_flag: This flag controls whether to perform the backward sweep and use the DM-NRG 
method or not. Allowed tags are 'ON' or 'OFF'. In case of 'ON', the DM-NRG backward sweep 
is carried out and the complete set of eigenstates is used to calculate the spectral functions. In 
case of 'OFF' the usual NRG procedure is performed only. 

• text_swap_files_flag: This flag controls whether to generate various temporary scratch files 
as readable text files in addition to the binary files. Allowed tags are 'ON' or 'OFF'. 

• hoppings_on_site_energy_flag: This flag is important when energy-dependent density of 
states are considered. This flag controls whether the hoppings and the on-site energies are read 
from a file or not. If set to 'ON', the hoppings as well as the on-site energies must be computed 
using the 'he' utility that comes with the code. This utility generates some data files in the 
. /dos_mapping folder, which are next read by the fnrg binary. When a flat band is considered, 
this flag can be set to 'OFF' and then the hoppings are generated on the fly. 

4.3.3 <SECTION-SYMMETRIES> 

In this section the types of symmetries used can be set. As many symmetries must be specified 
as were set by the 'symmetry _no' variable. Each line must contain the name of a symmetry, 
and a minimum and maximum value for the corresponding representation index. Those states 
for which these bounds are exceeded are simply discarded, and to avoid such uncontrolled 
truncation, a large value (~ 20) should be used for most symmetries. The possible symmetry 
types and the corresponding tags that the present version of the code can handle are as follows: 
'£/(l)', , SU(2)\ '^(2)', 'charge-SU(2)\ Other symmetries will also be included in later versions 
of the code. 

4.3.4 <SECTION-BLOCK_STATES> 

In this section the quantum numbers of the block states are set. A matrix of dimensions 
block_state-.no x symmetry jno with integer entries enumerating the representation indices of the 
block states must be provided. The number of rows has to be equal to the value 'block_state_no' 
set in < SECTION- PARAMETERS >. Integers in a row u (u = 1, . . . , block, state .no) corre- 
spond to the representation indices (quantum numbers) of the state u. Thus each row contains 
symmetry .no integer entries. Note that it is necessary to provide the double of the usual quan- 
tum numbers in case of SU (2) symmetries in order to be able to work with integer numbers. 
It is important that the states must be ordered by symmetries: The ordering is done with 
respect to the first symmetry, then with respect to the second symmetry used and so on. While 
reading the input file, the code checks if the quantum numbers are in ascending order, and if 
not, then an error message is generated and the run stops. 

4.3.5 <SECTION-LOCAL_STATES> 

In this section the quantum numbers are set. An integer type matrix of dimensions lo- 
cal_state_noxsymmetry_no must be provided. The number of rows must be equal to the value 
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'local_state_no' set in < SECTION- PARAMETERS >. Otherwise, the rules are the same as for 
<SECTION-BLOCK_STATES>. 

4.3.6 <SECTION-LOCAL_STATES_SIGNS> 

In this section the signs of the local states must be set, and a vector of length 'local_state_no', 
containing ±l-s must be provided. The sign of the local state is -1 / +1 if it contains an odd 
/ even number of fermions. 

4.3.7 <SECTION-BLOCK_HAMILTONIAN> 

This section contains the definition of the interaction Hamiltonians, Ti$. 

• <BLOCK_HAMILTONIAN_TERM> tag denotes the beginning of the definition of 
a new interaction Hamiltonian term. 

• <BLOCK_HAMILTONIAN_NAME> is a character string of maximum length char [256] 
to define the name of the Tig term. It could be any string. 

• <BLOCK_HAMILTONIAN_COUPLING> a rcal*8 number, which sets the cou- 
pling strength, (G$), corresponding to Tig. In this line you also specify the name (nota- 
tion) of the coupling. By writing a line "J = 0.736", e.g., you teach the code that this 
coupling is called J, and its value is 0.736. 

• <BLOCK_HAMILTONIAN_REPRESENTATION_INDEX> 

the quantum numbers of the Hamiltonian term, an integer vector of length 'symmetry _no'. 
Under any circumstances the quantum numbers must be for all the symmetries, other- 
wise the problem is badly posed. The representation index of the Hamiltonian may look 
therefore redundant at first, but to have a uniform treatment of all operators in the input 
file, we decided to keep this line. 

• <BLOCK_HAMILTONIAN_MATRIX> This is a matrix of dimension block_state_no 
x block_state_no. This matrix contains the matrix elements of the Hamiltonian Hg. 
between the states (multiplets) enumerated in <SECTION-BLOCK_STATES>. Some 
Mathematica 6.0 notebooks are provided in directory mathematicajfiles which can be 
used to generate these matrices in a routine way. 

The above series of entries must be repeated for every local Hamiltonian term, H$. The number 
of these structures must be equal to coupling_no defined in < SECTION- PARAMETERS >. 

4.3.8 <SECTION-HOPPING_OPERATORS> 

This section provides information on the hopping operators. As a rule, we provide information 
only on creation operators. No information on annihilation operators is needed. 

• <HOPPING_OPERATOR> tag denotes the beginning of the definition of a new 
hopping operator. 

• <HOPPING_OPERATOR_NAME> is a character string with maximum length char[256] 
to define the name of the corresponding hopping operator. (Could be something like 
"spin_up_hopping" , e.g.) 

• <HOPPING_OPERATOR_REPRESENTATION_INDEX> the quantum numbers 
of the hopping operator, an integer vector of dimension 1 x 'symmetry _no' must be given. 
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• <HOPPING_OPERATOR_SIGN> The sign is -1 for a fermionic operator and 1 for 
a bosonic operator. 

• <HOPPING_OPERATOR_MATRIX> Contains the reduced matrix elements of the 
hopping operator (~ / f ) between the block states specified in <SECTION-BLOCK_STATES>. 
It is a block_statc_no x block_statc_no matrix. Wc provided a few Mathematica note- 
books which you can use to generate these matrix elements in a simple way. 

The above series of entries must be repeated as many times as many hopping operators one 
has. The number of these structures must thus agree with the 'hopping_operator_no' set in 
<SECTION-PARAMETERS>. 

4.3.9 <SECTION-SPECTRAL_OPERATORS> 

This section provides information on the spectral operators, i.e., operators, whose spectral 
functions are computed. As a rule, we provide information only on the 'creation operators', 
and no information on annihilation operators is given (see Section [2.21 for a precise definition). 
The following structure must be repeated for every spectral operator. Clearly, the number of 
structures must be equal to the parameter spectral_operator_no set in section < SECTION- 
PARAMETERS >. 

• <SPECTRAL_OPERATOR> tag denotes the beginning of the definition of a new 
spectral operator. 

• <SPECTRAL_OPERATOR_NAME> is a character string with maximum length 
char[256] to define the name of the corresponding spectral operator. (E.g., "spin opera- 
tor") 

• <SPECTRAL_OPERATOR_REPRESENTATIONJNDEX> the quantum num- 
bers of the spectral operator. An integer like vector of dimension 1 x symmetry _no is 
expected. 

• <SPECTRAL_OPERATOR_SIGN> The sign is -1 for a fcrmionic-likc operator and 
1 for a bosonic operator. 

• <SPECTRAL_OPERATOR_MATRIX> A matrix of dimensions blockj3tatc_no x 
block_state_no is expected, containing the reduced matrix elements of the spectral operator 
between the block states specified in <SECTION-BLOCKJ3TATES>. We included in this 
package a few examples of Mathematica notebooks, which can be used to generate these 
matrices relatively easily. 

4.3.10 <SECTION-STATIC_OPERATORS> 

This section provides information on the operators for which static averages are evaluated. The 
following structure must be repeated for every static operator. 

• <STATIC_OPERATOR> tag denotes the beginning of the definition of a static oper- 
ator. 

• <STATIC_OPERATOR_NAME> is a character string with maximum length char[256] , 
defines the name of the static operator being defined. 

• <STATIC_OPERATOR_REPRESENTATIONJNDEX> quantum numbers of the 
static operator, an integer type vector of dimension 1 x symmetry _no is expected. 
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• <STATIC_OPERATOR_SIGN> This sign is -1 for a fermionic and 1 for a bosonic 
operator. 

• <STATIC_OPERATOR_MATRIX> Contains the reduced matrix elements of the 
static operator between the block states specified in <SECTION-BLOCK_STATES>. It 
is a block_state_no x block_state_no matrix. For the calculation of these matrix elements, 
one can use the mathematica codes provided with the package. 

The number of structures shown above must agree with the variable 'static_operator_no' set in 
<SECTION-PARAMETERS>. 

4.3.11 <SECTION-LOCAL_HAMILTONIAN> 

• <LOCAL_HAMILTONIAN_TERM> tag denotes the beginning of the definition of 
a local Hamiltonian, i.e. a term that occurs in Hi, 

local -coupling _no 

H. = 9a,i Hi 

a=l 

with Hi. a = Ti. a ({C_i \})- In other words, the local Hamiltonian Hi :Ct depends on the site 
index i only through the creation operators acting at site i, otherwise its structure is fixed. 
A trivial example of such operators are the energies, £a <-> g, <-» Hi tCt . In the 

present version of the code, we only handle uniform local couplings, g a ^ = g a , however, 
later versions shall also be able to treat a user-specified coupling set, g a ^. Note that the on- 
site energy terms are specified under the section <SECTION-LOCAL_HOPPING_OPERATORS>, 
not here. Here you can define additional local terms in the Hamiltonian, e.g., supercon- 
ducting pairing terms. 

• <LOCAL_HAMILTONIAN_NAME> is a character string with maximum length 
char [256], defines the name of the Hi term. 

• <LOCAL_HAMILTONIAN_COUPLING> is a rcal*8 number to set the coupling 
strength, g a , corresponding to the local term, Hi ta - In the present version of the code the 
coupling strength does not change along the Wilson chain. 

• <LOCAL_HAMILTONIAN_REPRESENTATIONJ[NDEX> the quantum num- 
bers of the Hamiltonian term, an integer type vector with dimension 1 x 'symmetry _no'. 
Similar to the Hamiltonian, these quantum numbers must be all 0, otherwise the Hamil- 
tonian is not invariant under the symmetry assumed. 

• <LOCAL_HAMILTONIAN_MATRIX> It is a localj3tate_no x local_state_no ma- 
trix, defined on the basis given in <SECTION-LOCAL_STATES>. This contains the 
(reduced) matrix elements of the coupling Hamiltonians, H a - 

4.3.12 <SECTION-LOCAL_HOPPING_OPERATORS> 

Here we specify matrix elements of the local hopping operators between the local states enu- 
merated in <SECTION-LOCAL_STATES>. As a rule, we provide the information only for the 
creation operators. 

• <LOCAL_HOPPING_OPERATOR> denotes the beginning of the definition of a 
hopping operator. 

• <LOCAL_HOPPING_OPERATOR_NAME> is a character string with maximum 
length char[256]. Defines the name of the local hopping operator. 
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• <LOCAL_HOPPING_OPERATOR_REPRESENTATION_INDEX> the quan- 
turn numbers of the local hopping operator, an integer type vector with dimension 1 x 
'symmetry _no' must be given. 

• <LOCAL_HOPPING_OPERATOR_SIGN> This sign is -1 for a fermionic and 1 for 
a bosonic operator. 

• <LOCAL_HOPPING_OPERATOR_MATRIX> Contains the reduced matrix ele- 
ments of the static operator between the local states specified in <SECTION-LOCAL_STATES>. 
It is a local_state_no x local_state_no matrix. For the calculation of these matrix elements, 

one can use the mathematica codes provided with the package. 

• <LOCAL_ON_SITE_ENERGY_MATRIX> Contains the reduced matrix elements 

of the on-site energy term between the local states specified in <SECTION-LOCAL_STATES>. 
It is a local_state_no x local_state_no matrix. For the calculation of these matrix elements, 
one can use the mathematica codes provided with the package. 

The number of structures shown above must agree with the variable 'hopping_opcrator_no' set 
in the <SECTION-PARAMETERS>. 

4.3.13 <SECTION-SPECTRAL_FUNCTION> 

In this section one specifies the pairs of spectral operators, whose spectral function is computed. 
This section must contain 'spectral_function_no' rows, each of them corresponding to a given 
pair of operators. 

To compute the spectral function of the Green's function of the form 

G% Bi (t) = -i([A(t),B^(0)} i )Q(t), (4.3) 
one must add a row of the form 

{name_of_^; name_of_-B' } , 

where 'name_of_A^' and 'name. of_Bt> are the names of the 'creation operators' set up under 
<SECTION-SPECTRAL_OPERATORS>. 

4.3.14 <SECTION-SPECTRAL_FUNCTION_BROADENING> 

Here you specify the parameters of the broadening, used to evaluate the spectral functions. 
Note that the code fnrg does not generate the spectral functions, you need to run the utility 
sfb for that immediately after the run. This part of the input file is thus used by sfb. Note 
that one can improve/change the broadening parameters after running the code too, and then 
run sfb without running fnrg again. For details on the broadening, see Section f2~2|) . 
The following parameters must be set in this section: 

• broadening_method: the method that is used for the broadening of the spectral func- 
tion along the real axis. Two methods are available: LOG_GAUSS and INTERPOLA- 
TIVE.LOG_GAUSS. 

• broadening_parameter: The value of the broadening parameter (half-width of the log 
Gauss function in both the log-Gauss and interpolative log-Gauss type of broadening). 
This parameter is used no matter what method is used. 

• broadening_energy_minim: The minimum energy for which the broadening is per- 
formed (broadening _energyjminim < \uj\ < broadening jenergy -maxim). This parame- 
ter is used for both methods. This quantity must be positive. 
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broadening_energy_maxim: The maximum energy for which the broadening is per- 
formed (broadening -energy ..minim < \u>\ < broadening jenergyjmaxim). This quantity 
must be positive. 

broadening_grid_mesh: The number of points along the real axis, for which the broad- 
ened spectral function is computed. The binary sfb generates a logarithmic mesh between 
[-broadening_energy_maxim, -broadening_energy_minim] and [broadening_energy_minim, 
broadening_energy_maxim] for which the spectral function is computed. 

quantum_temperature: This parameter is used only with the INTERPOLATIVE_LOG_GAUSS 
method. It sets the scale below which a linear interpolation is used rather than a loga- 
rithmic interpolation. 

spectral_function_nrg_even: This is a flag that is relevant for NRG runs only. Allowed 
tags are 'YES' or 'NO': In case of 'YES' the spectral function is computed from the even 
NRG iterations, and a corresponding data file is generated. If set to 'NO', the code skips 
this calculation. 

spectral_function_nrg_odd: This is a flag that is relevant for NRG runs only. Allowed 
tags are YES' or 'NO': In case of YES' the spectral function is also computed from odd 
NRG iterations. If set to 'NO', the code skips this calculation. 

spectral_function_dmnrg: This is a flag parameter that controls whether the spectral 
function from the full NRG spectrum using the density matrix formalism is computed. 
Allowed tags are YES' or 'NO'. In case of YES' the DM-NRG calculation is also per- 
formed. 

static_average_dmnrg: Flag parameter that controls whether the static average of the 
form (A^ B) is also computed by sfb for the two spectral operators, in addition to the 
calculation of the spectral functions (see the notations in the previous section). Allowed 
tags are YES' or 'NO'. 

green_function_nrg_even: This flag controls whether the real part of the Green's func- 
tion is computed from the even part of the spectral function generated from the nrg data, 
by performing a Hilbcrt transform. Allowed tags are YES' or 'NO'. The real part of the 
Green's function can be computed only after the calculation of the spectral function, i.e. 
if the spectral_function_nrg_even flag is set to YES. 

green_function_nrg_odd: This flag controls whether the real part of the Green's func- 
tion is computed from the odd part of the spectral function generated from the nrg data, 
by performing a Hilbcrt transform. Allowed tags are YES' or 'NO'. The real part of the 
Green's function can be computed only after the calculation of the spectral function, i.e. 
if the spectral_function_nrg_odd flag is set to YES. 

green_function_dmnrg: This flag controls whether the real part of the Green's function 
is computed from the dmnrg spectral function by performing a Hilbcrt transform. Allowed 
tags are YES' or 'NO'. The real part of the Green's function can be computed only after 
the calculation of the spectral function, i.e. if the spectral_function_dmnrg flag is set 
to YES. 

green_function_grid_mesh: The number of points on a log-mesh used for doing the 
Hilbcrt transform. Once the spectral function for a combination of spectral opera- 
tors is computed, the real part of corresponding Green's function is generated by per- 
forming a Hilbcrt transformation. To compute the Hilbcrt transform, a logarithmic 
mesh is generated. The energy limits for the real/imaginary part of the Green's func- 
tions are the same as those for the spectral function , i.e. [-broadening_energy_maxim, 
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-broadcning_cnergy_minim] and [broadcning_cncrgy_minim, broadcning_energy_maxim] . 
The energy values for which the Green's function is computed correspond to those for 
the spectral function. Cubic splines are used for the interpolation of the spectral func- 
tion. Therefore, for too small green_function_grid_mcsh values sometimes oscillations are 
observed in the energy dependence of the real part of the Green's function. In this case, 
green-function_gridjmesh must be increased until the oscillations vanish. Increasing this 
parameter has only a small impact on computing time. 
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4.4 Outputs of the DM-NRG code 

4.4.1 Output directory structure 

The output of the code has the following directory structure: 



. /results/ 
. /results/ 



automat ic_directory_name/ 
(code generated) 
log_f ile . dat 



input . dat 

spectral_function_dmnrg_* 
spectral_function_even_* 
spectral _function_odd_* 
green_f unct ion_dmnrg_* 
gr een_f unct ion_e ven_* 
gr een_f unct ion_odd_* 
static_average_dmnrg* 
energy_levels_even* 
energy_levels_odd* 
data/ 



all results of the code 
all data for a given run 



spectral_operators/ spectral_operator_* . dat 
static_operators/ static_operator_* . dat 
sectors/ sectors.* . dat 
transformations/transformation.* . dat 
states/ states.* . dat 
delta_peaks/delta_peaks_* .dat 
spectral_f unctions/ spectral_f unction_* . dat 
mapping/dos .dat 
mapping/hopping_couplings* . dat 
mapping/on_site_energies* .dat 
mapping/density_of .states* .dat 
m_star.dat 



Table 4.1: Output directory structure of the code. Directory automatic_directory_ name/ is 
generated automatically from the parameters in the input file. The "*" in the file names also 
stands for a list of the parameters (names of the couplings etc.) generated automatically. 

All outputs of the code, which are saved on disk, are located in the directory results. In 
this directory the program generates automatically a directory corresponding to the parame- 
ters of the given run. The name of this directory, automatic_directory_name/, contains the 
values of the following input variables: model, lambda, maxstatejrio, iteration jno , temperature, 
symmetry_.no , symmetry J,ypes, coupling_no, and a file log_file.dat containing the log/error 
messages during the run. 

The automatic_directory_name/, folder has the following structure: 

• input.dat is a copy of the original input file. 

• Various *.dat files contain the final results for the spectral functions, static averages, 
green's functions and the energy levels. The names of these files are descriptive, and 
contain some additional tags such as, e.g., the names of quantum operators for which the 
spectral functions have been computed. 

After running fnrg, only a directory data containing the raw data is generated. The spectral 
functions, Green's functions, static averages files are generated only later by running the sfb 
utility. The energy level files are also generated separately, by running the es utility. 

• Directory data contains the data files generated by fnrg. Since there are hundreds of files, 
this directory is divided into subdirectories with descriptive names and containing the following 
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files: 

mapping/ 

sectors/ sectors.* . dat 

states/ states.* . dat 

delta_peaks/delta_peaks_dmnrg_* 

delta_peaks/delta_peaks_nrg_* 

spectral_operators/ spectral_operator_* 

static_operators/static^average_dmnrg_* 

static_operators/static^average_nrg_even_* 

static_operators/static^average_arg_odd_* 

transformations/transformation.* . dat 

spectral Junctions/ spectral Junction.* . dat 

m_star.dat Here again, *'s stand for operator names and code-generated tags such as 
the iteration number etc. 

4.4.2 Detailed description of the output files in the folder 

. / results/ automat ic_directory_name/ data 

This folder contains the original data generated/used by the code fnrg and the utilities he and 
sfb. It contains the following files/folders: 

• mapping/ : This folder contains information on the hopping parameters and on site energies, 
in case a user-defined density of states has been used. This directory is a copy of the directory 
dos_mapping, located in the main directory. It contains the following files: 

• mapping/hopping_couplings*.dat: Contains the data for the hopping couplings along 
the chain as computed by the he (hopping energies) utility. Each line contains information for 
one iteration, the first column being the iteration number and subsequent columns representing 
the values of the hopping amplitudes along the Wilson chain. Each column corresponds to a 
hopping operator, so the number of columns has to agree with the 'hopping_operator_no'. The 
tag "self-consistency" in the file name has no relevance here, it was included for later purposes. 

• mapping/on_site_energies*.dat contains the data for the on-site energies along the chain 
as computed by the he (hopping energies) utility. Each line contains information for each 
iteration, the first column being the iteration number and the subsequent columns being the 
values of the on-site energies. Each column corresponds to a hopping operator, so the number 
of columns has to agree with the 'hopping_operator_no'. The tag "self-consistency" in the file 
name has no relevance here, it was included for later purposes. 

• sectors/sectors_(iteration_no).dat The code is organized using so-called sectors. A given 
sector is characterized by a list of representation labels, specifying the symmetries of the states 
of that sector. Similarly, operators are also organized using the notion of sectors. The files 
sectors_(iteration_no) .dat contain the information on the size and symmetry of sectors in 
iteration iteration -no. 

This file has two separate parts. The first part describes the sectors in the new basis at iteration 
iteration_no. The first row of the file contains the total number of the sectors in the new basis 
(new -blocks ectorjno) . The subsequent new-blocksector-no lines contain the information for 
each sector (from to new-blocksector-no-1): the representation indices of the sectors, the 
dimensions of the sectors, the first states of the sectors among the list of all the basis states 
(multiplets), and the degeneracies of the sectors. Note that the sum of the dimensions multiplied 
with the corresponding degeneracy values must be equal to the total number of basis states of 
the new block state basis. 
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The second part of the files sectors_(iteration_no) .dat contains the same information on 
the states kept from the previous iteration (iteration iterationjno- 1) . 

• states/states_(iteration_no).dat: These files contain the information for the new basis 
states at iteration iterationjno, and the states kept from the previous iteration (iteration iter- 
ation_no-l). States know about their parent states, and this information is also displayed in 
these files. 

The first row of the file contains the total number of the new basis states (newJAockstatejno) 
The next newJAockstatejno lines contain the information on the multiplets after the diagonal- 
ization has been done. Each line contains the number (label) of the sector the state belongs 
to, the index of the state within this sector, the list of representation indices, the labels of 
the parent and local states the state was constructed from, the sign of the state (this is for 
block states), whether the state is kept (K) or eliminated (E) after truncation, and finally the 
eigen-energy of the state. Note that the sector number cannot be less than zero and it must be 
smaller than newJ)lock_sector-no, specified in file sectors-(iteration_no) .dat. Also the variable 
sector Jmde.x must be non-negative and less than the dimension of the sector given also in file 
sectors_(iteration_no) .dat. 

The second part of the file states/states_(iteration_no) .dat gives the same information 
on the states kept from the previous iteration. 

• transformations/transformation_(iteration_no).dat: These files contain the informa- 
tion on the matrix O n that is used to diagonalize the Hamiltonian, and to transform the 
hopping operators at site iteration-no, and the spectral operators (acting at site 0) to the new 
basis as 

A n -> O n A n Ol (4.4) 

These files are only generated if the flag text_swap_files_flag in the file input . dat is set ON. 
The first row of the file contains the total number of the sectors in the new basis in iteration 
iterationjno, and this is followed by three terms for each sector: the sector number, the size of 
the sector, and finally the real*8 matrix of size sector_size x sector_size, transforming the given 
sector. 

• transformations/transformation_binary_(iteration_no).dat contains the information 
for the O matrix performing the diagonalization at iteration iterationjno in a binary form. 
These files are only generated if the flag binary_swap_files_flag in the file input . dat is set to 
ON. 

• spectral_operators/spectral_operator_(operafor_name)_(iteration_no).dat contains 
the matrix elements of a spectral operator acting at site of the chain (for example d Q or 
S z ) in a sector-decomposed form, after the diagonalization in iteration step iterationjno has 
been done. These operators (which are 'creation operators') are defined and initialized in the 
input.dat file, <SECTION-SPECTRAL_OPERATORS>. These files are only generated if 
the flag text_swap_files_flag in the file input . dat is set ON. 

The matrix of a spectral operator can be decomposed to sector_no by sector _no small blocks. 
Depending on the quantum numbers of the operator, several sectors (actually most of them) 
contain only zero elements. These are considered as non-active sectors. The first line of the file 
contains the total number of the active sectors (those with proper change of quantum numbers). 
Then for each active sector the row and column index of the sector is given and is followed by 
the row and column dimension of the sector, and the matrix elements of the operator in that 
sector. 

• spectral_operators/spectral_operator_binary_( operator name)_(iteration_no).dat: 

These files contain the same information as the files spectral_operator_( operat or name) j*. dat, 
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excepting that the data are in a binary form. These files are only generated if the flag bi- 
nary_swap_files_flag in the file input . dat is set to ON. 

• static_operators/static_operator_( operator name)_(iteration_no).dat: They contain 
the matrix elements of a static operator acting at site of the chain (for example no) in a 
sector-decomposed form, after the diagonalization in iteration step iteration_no has been done. 
These operators (which are 'creation operators') are read from the input . dat file, <SECTION- 
STATIC_OPERATORS>. These files arc only generated if the flag text_.swap_files.flag in the 
file input . dat is set ON. 

Static operators must have a block-diagonal structure by symmetry, therefore only these sectors 
are active. These files are organized otherwise in the same way as the files spectral_operator_* . 

These files are only generated if the flag text_swap_files_flag in the file input . dat is set to ON. 
Furthermore corresponding binary files are also generated if the flag binary_swap_file.s_flag in 
the file input . dat is set to ON. 

• delta_peaks/delta_peaks_dmnrg_( operator l)_(operator 2) _(iteration_no ).dat con- 
tain the raw spectral function peaks for the spectral function of operator 1 and operator 2, 
generated from iteration iteration_no. These files contain two columns: the first one is the 
energy and the second one the contribution to the amplitude of the spectral function at that 
energy in iteration iteration_no. The calculation is performed at a temperature T specified 
in the input file, input.dat. These output files are used by the utility sfb to compute the 
DM-NRG spectral functions, combining the data from all iterations. 

• delta_peaks/delta_peaks_nrg_( operator name l)_(operator name 2) _(iteration_no) 

.dat contains peaks for the spectral function of operator 1 and operator 2, generated from 
iteration iteration_no. These files contain two columns, the first one being the renormalization 
energy and the second one the amplitude of the spectral function at the corresponding energy 
and iteration iteration_no. The calculation is performed at a temperature T specified in the 
input file, input.dat. These files are used to compute the NRG spectral function combining 
data from either even or odd iterations. 

• spectral_functions/spectral_function_dmnrg_( operator l)_(operator 2) 
_(broadening_par) _(iteration) .dat contains two columns, the first one being the energy and 
the second one the real axis value of the spectral function corresponding to the two operators 
from the filename at the iteration iteration_no. Data for negative as well as positive energies are 
saved into this file. These particular files are relatively important for analyzing the contribution 
of each iteration to the total spectral function. 

• static_operators / static_average_nrg_even_(operator_name) .dat and 
static_average_nrg_odd_(operator_name).dat: These files contain average of a static op- 
erator as a function of temperature, computed by the standard NRG prescription, from the 
even/odd iterations, respectively. 

• static_operators/static_average_dmnrg_(operator_name).dat and 
static_average_dmnrg_(operatorl)_(operator2).dat: These files contain the expectation 
values of the static operators and of the pairs of spectral operators specified in the file input . dat, 
computed using the DM-NRG procedure at a temperature T, also specified in the file input . dat. 

• m_star.dat: m_star is the value of the iteration_index at which truncation of the states begins 
(new_state_no > max_state_no) , i.e., for iterations larger then m_star the states, transformation 
matrices, spectral operators have to be saved. m_star remains fixed for the rest of the NRG 
steps. This file is set on input and by function get_m_star () in cut_of f . cc and read by various 
functions when spectral functions are generated later. 
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4.4.3 Analyzing the data 



There are several utilities that we provide to analyze the data you generated. 

• First of all, spectral functions and their Hilbert transforms can be generated/analyzed by 
the utility sfb. The broadening parameters must be set in the file input . dat sitting in 
the main folder. You can set there the broadening method, as described in Chapter 14.31 
(See also the last section of Chapter [2]) Depending on the parameters set this utility 
generates the following files. 

spectral_function_dmnrg_( operator 1) -(operator 2) -( broadening), dat 

These files contain the broadened spectral function of operators operatorl and opera- 
tor. Correspondingly, in case the Hilbert transform flag is ON in the file input . dat, 

the real part of the corresponding Green's function is also generated in file green_function_dmnrg_*.dat. 

spectral_function_even_( operatorl) _(operator2) -(broadening) . dat and 
spectral_function_odd_( operatorl) -(operator 2) -(broadening) . dat: 

These files contain the broadened spectral functions of operatorl and operator2 7 as 
computed by the traditional NRG method. Data arc collected from even and odd 
iterations respectively If the Hilbert transform flag is ON in the file input.dat, 
then the real parts of the corresponding Green's functions are also generated and 
stored in the files green_function_even_*.dat and green_function_odd_*.dat. 

• The energy spectrum (finite size spectrum) can be analyzed by running the utility es. 
While running this utility you can specify quantum number filtering as well as the num- 
ber of states that you want to keep track of. This is a useful utility to obtain the fixed 
point spectrum of a model, identify the corresponding conformal field theory, the di- 
mensions of the irrelevant operators, or to extract Fermi liquid parameters. The results 
are stored separately for even and odd iterations in the files energy _levels_even* and 
energy _levels_odd*. 
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Appendix A 



Input file for the Kondo model 

with ^charge! 1 ) x ^spin(l) 

symmetries 



# The comments must start with the symbol # . 

# Empty lines are ignored. 

######################################################### 
< SECTI ON-P ARAMETERS > 

######################################################### 

# Explanation of variables 

# 'model' represents the type of the model that is 

# used for the calculations. It can be any name. 

# 'lambda' represents the value for the logarithmic 

# discretization parameter. Usually it is fixed between 2 and 3. 

# 'max_state_no ' is the number of kept states 

# after an iteration. 

# ' iteration_no ' is the number of iterations performed 

# along the NRG run. Typically it is set to 40 - 70. 

# ' symmetry_no ' specifies how many symmetries 

# are used. The symmetries are characterized 

# in <SECTION-SYMMETRIES> 

# ' coupling_no ' specifies how many couplings 

# are present in the interaction Hamiltonian. 

# 'local_coupling_no' is the number of 

# Hamiltonians at the sites we add to the Wilson chain 

# 'spectral_operator_no' is the number of spectral operators 

# acting on site 0, we keep track of. 

# These spectral operators are characterized 
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# in section <SECTION-SPECTRAL-OPERATORS> 



# ' static_operator_no ' is the number of static operators 

# acting on site 0, we keep track of. 

# These static operators are characterized 

# in section <SECTI0N-STATIC-0PERAT0RS> 

# 'hopping_operator_no ' is the number of hopping operators. 

# These hopping operators are characterized 

# in section <SECTI0N-H0PPING-0PERAT0RS> 

# 'block_state_no ' specifies the number 

# of initial states/multiplets . These states are characterized 

# by the representation indices and are listed 

# in section <SECTI0N-BL0CK_STATES> 

######################################################### 
######################################################### 

model = kondojnodel # the name of the present model 
lambda = 2.0 # discretization parameter 
max_state_no = 100 # maximum number of kept multiplets 
iteration_no = 20 # allowed number of iterations 
symmetry_no = 2 # the number of symmetries 

coupling_no = 2 # no of Hamiltonian terms, i.e. J and B in this case 

spectral_operator_no = 3 # no of operators for which correlation functions are computed 
static_operator_no = 1 # no of operators for which the static average is computed. 
hopping_operator_no = 2 # the number of hopping operators 
block_state_no = 8 # the number of initial block states 
local_state_no = 4 # the number of local states added in each iteration 
local_coupling_no = 1 # the number of Hamiltonian terms on the added site 
spectral_function_no = 3 # the number of spectral functions that need to be generated 
interval_no = 1000 # mash grid between the maximum and minimum energy at each iteration 
degeneracy_threshhold = le-6 # energy thresh-hold for discarding states 
temperature = 0.0 # temperature used for the run 

######################################################## 
< / SECTI 0N-P ARAMETERS > 
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######################################################## 



######################################################## 

<SECTION-FLAGS> 

######################################################## 



# dmnrg_flag controls whether the backward procedure 

# has to be done and the full set of eigen-states be 

# used for the calculation of the spectral function. 

# ON - backwards iteration is done. 

# OFF - only the nrg calculation is performed. 

# text_swap_f iles_f lag = flag that fixes if the text mode is 

# used to save the files. The files will be saved in binary mode anyway. 

# ON - the files are saved in the text mode and kept on the disk. 

# 0FF= the files are not saved in the text mode. 

# binary _swap_f iles_f lag = flag that controls whether the 

# unnecessary binary files will be removed or not not after the calculation is done. 

# ON - the files will be kept on the disk. 

# OFF- the files will be removed from the disk. 

# hoppings_on_site_energies_f lag = flag that controls 

# the reading of the hoppings and the on-site energies. 

# ON - the hoppings and the on site energies 

# are read from the files in the results/mapping/ folder 

# These files are generated by the he utility. 

# OFF - the hoppings are computed on the fly 

# assuming a flat density of states (DOS =0.5) on [-1, 1], 

# and the on-site energies are set to zero. 

dmnrg_flag = ON 
text_swap_f iles_f lag = OFF 
binary_swap_f iles_f lag = ON 
hoppings_on_site_energies_f lag = OFF 

######################################################## 

< / SECTION-FLAGS > 

######################################################## 



######################################################### 

<SECTION-SYMMETRIES> 

######################################################### 

# As many symmetries must be defined below 

# as set by the variable ' symmetry_no ' 

# in the previous section. 

# The possible symmetry types are the following: 
#U(1), SU(2), Z(2), charge_SU(2) . 

# 

# Next the limits for the representation index 
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# are set . All states with smaller/larger 

# representation indices will be cut off 

# For U(l) symmetries the representations are labeled 

# by integers (2 S_z) and for SU(2) by non-negative 

# integers (2S) . 

# For the Z(2) symmetry the representations are labeled 

# by integers 

# - even parity states 

# 1 - odd parity states 

######################################################### 



# charge symmetry corresponding to Q 
U(l) -10 10 

# spin symmetry corresponding to 2 S_z 
U(l) -20 20 



######################################################### 

</SECTION-SYMMETRIES> 

######################################################### 



######################################################### 

<SECTI0N-BL0CK_STATES> 

######################################################### 

# This section gives the representation indices of the 

# block states of the first iteration. 

# The number of rows has to be equal 

# with the 'block_state_no ' in the first section, 

# and the number of columns to the ' symmetry_no ' . 

# Each column contains indices corresponding to a symmetry. 
######################################################### 

-1 -1 
-1 1 
-2 



2 

1 -1 
1 1 

######################################################### 
</SECTION-BLQCK_STATES> 

######################################################### 



######################################################### 
< SECTI 0N-L0CAL_STATES > 

######################################################### 

# This section gives the representation indices 

# of the local states . The number of rows has to be equal 

# with the ' local_state_no ' in the first section 
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# and the number of columns to the ' symmetry_no ' . 

# Each column contains indices corresponding to a symmetry. 
######################################################## 

-1 
-1 

1 

1 

######################################################## 
</SECTION-LOCAL_STATES> 

######################################################## 



######################################################## 
<SECTION-LOCAL_STATES_SIGNS> 

######################################################## 

# These are the signs of the added states. 

# The number of signs has to be equal with 

# the number of ' local_state_no ' . 

# All the signs must be enumerated in one line . 

# For even electron number, the sign is 1, 

# while for odd electron number the sign is -1. 
######################################################### 

1-1-11 

######################################################### 
< /SECTION-LOCAL_STATES_SIGNS > 

######################################################### 



######################################################### 
<SECTIQN-BLOCK_HAMILTONIAN> 

# The total Hamiltonian will be written as a sum 

# of the Hamiltonian terms. 

tlo — Os-tls 

# There has to be as many <HAMILTONIAN_TERM> items 

# as 'coupling_no' were set in the first section. 

# The first tag is the name of the Hamiltonian term, 

# the second one gives the representation indices and 

# at last comes the matrix elements that will be 

# read into the sectors. 



# first term in the Hamiltonian 

<BLDCK_HAMILTONIAN_TERM> 

<BLQCK_HAMILTQNIAN_NAME> 
HJCondo 

</BLOCK_HAMILTQNIAN_NAME> 
<BLOCK_HAMILTONIAN_CQUPLING> 

# J - exchange coupling 
J = 0.5 



■ r jl 



< /BLOCK_HAMILTONI AN_COUPLING> 
<BLQCK_HAMILTQNIAN_REP RESENT ATION_INDEX> 


</BLOCK_HAMILTONIAN_REPRESENTATION_INDEX> 



<BLQCK_HAMILTQNIAN_MATRIX> 

























































0.25 


























-0.75 


























0.25 


























0.25 

























































</BLOCK_HAMILTONIAN_MATRIX> 



< /BL0CK_HAMILT0NI AN_TERM> 

# second term in the Hamiltonian 

<BL0CK_HAMILT0NIAN_TERM> 

<BL0CK_HAMILT0NIAN_NAME> 
H_Zeeman 

</BLOCK_HAMILTQNIAN_NAME> 
<BL0CK_HAMILT0NIAN_C0UPLING> 

# B - external magnetic field 
B = 0.01 

</BL0CK_HAMILT0NIAN_C0UPLING> 
<BL0CK_HAMILT0NIAN_REP RESENT ATI0N_INDEX> 


</BLOCK_HAMILTONIAN_REPRESENTATION_INDEX> 



<BLOCK_HAMILTONIAN_MATRIX> 



-0.5 


























0.5 


























-0.5 





























0.5 




















0.5 





























0.5 


























-0.5 


























0.5 



</BLOCK_HAMILTONIAN_MATRIX> 
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</BLOCK_HAMILTONIAN_TERM> 



######################################################## 

</SECTION-BLOCK_HAMILTONIAN> 

######################################################## 
######################################################## 

<SECTION-HOPPING_OPERATORS> 

# There must be as many <HOPPING_OPERATOR> items 

# as 'hopping_operator_no ' were set in the first section 

# The matrices must contain the Reduced Matrix Elements. 
######################################################## 

# here is the first hopping operator 
<H0PPING_0PERAT0R> 
<HOPPING_OPERATOR_NAME> 

f _N_dagger_up 

</HOPPING_OPERATOR_NAME> 

<HOPPING_OPERATOR_REPRESENTATION_INDEX> 



</H0PPING_0PERAT0R_REP RESENT ATION_INDEX> 
<HOPPING_OPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 



1 1 



-1 



</HOPPING_OPERATOR_SIGN> 



<HOPPING_QPERATOR_MATRIX> 















1 













1 


















0.707106781186 0.707106781186 

























-0.707106781186 
0.707106781186 









</H0PPING_0PERAT0R_MATRIX> 



</H0PPING_0PERATQR> 
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# here is the second hopping operator 
<HOPPING_QPERATOR> 
<HOPPING_QPERATOR_NAME> 

f _N_dagger _down 
</HOPPING_OPERATOR_NAME> 
<HOPPING_QPERATOR_REPRESENTATION_INDEX> 
1 -1 

</HQPPING_OPERATOR_REP RESENT ATION_INDEX> 
<HOPPING_QPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</HOPPING_OPERATOR_SIGN> 



<HOPPING_OPERATOR_MATRIX> 



















































1 


























0.707106781186 























0.707106781186 





















































0.707106781186 


-0.707106781186 


























-1 









</H0PPING_0PERAT0R_MATRIX> 



</H0PPING_0PERATQR> 



######################################################## 
</SECTI0N-H0PPING_0PERAT0RS> 

######################################################## 



######################################################## 
<SECTI0N-SPECTRAL_0PERAT0RS> 

######################################################## 

# In this section we set up the matrices for the 

# operators for which the spectral function is computed 

# Matrices must contain the reduced matrix elements . 
####################################################### 

# here comes the first spectral operator. 
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< SPECTRAL_OPERATOR> 



<SPECTRAL_OPERATOR_NAME> 
S_z 

</SPECTRAL_OPERATOR_NAME> 

<SPECTRAL_OPERATOR_REP RESENT ATION_INDEX> 


</SPECTRAL_OPERATOR_REPRESENTATION_INDEX> 
<SPECTRAL_OPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
1 

< / SPECTRAL_OPERATOR_SIGN> 



< SPECTRAL_OPERATOR_MATRIX> 



-0.5 


























0.5 


























-0.5 





























0.5 




















0.5 





























0.5 


























-0.5 


























0.5 



< / SPECTRAL_0PERAT0R_MATRI X> 

< / SPECTRAL_0PERAT0R> 

# here comes the second spectral operator. 
<SPECTRAL_0PERAT0R> 

<SPECTRAL_0PERAT0R_NAME> 

f _0 .dagger _up 

</SPECTRAL_0PERAT0R_NAME> 

<SPECTRAL_0PERAT0R_REP RESENT ATI0N_INDEX> 
1 1 

</SPECTRAL_0PERAT0R_REPRESENTATI0N_INDEX> 
<SPECTRAL_0PERAT0R_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</SPECTRAL_0PERAT0R_SIGN> 

< SPECTRAL_0PERAT0R_MATRIX > 
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-0.707106781186 


() 




















0.707106781186 


























1 


























1 


























0.707106781186 


0.707106781186 












</SPECTRAL_OPERATOR_MATRIX> 



< / SPECTRAL_0PERAT0R> 

# here comes the third spectral operator. 

<SPECTRAL_0PERAT0R> 
<SPECTRAL_0PERAT0R_NAME> 
f _0 .dagger _down 
</SPECTRAL_0PERAT0R_NAME> 
<SPECTRAL_0PERAT0R_REP RESENT ATI0N_INDEX> 
1 -1 

</SPECTRAL_OPERATOR_REPRESENTATION_INDEX> 
<SPECTRAL_0PERAT0R_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operat 
-1 

< / SPECTRAL_0PERAT0R_SIGN> 



< SPECTRAL_0PERAT0R_MATRIX > 



















































1 


























0.707106781186 























0.707106781186 





















































0.707106781186 


-0.707106781186 


























-1 









< / SPECTRAL_0PERAT0R_MATRIX> 



< / SPECTRAL_0PERAT0R> 

######################################################## 
</SECTI0N-SPECTRAL_0PERAT0RS> 

######################################################## 



######################################################## 
<SECTI0N-STATIC_0PERAT0RS> 



5G 



######################################################## 

# In this section we set up the matrices for the 

# operators for which the static values are computed 

# These operators need to have representation 

# indices always zero ! 

####################################################### 

<STATIC_OPERATOR> 

<STATIC_OPERATOR_NAME> 

S_z 

</STATIC_OPERATOR_NAME> 

<STATIC_DPERATOR_REP RESENT ATION_INDEX> 


</STATIC_QPERATOR_REPRESENTATION_INDEX> 
<STATIC_OPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
1 

< /STATIC_OPERATOR_SIGN> 



<STATIC_OPERATOR_MATRIX> 



-0.5 


























0.5 


























-0.5 





























0.5 




















0.5 





























0.5 


























-0.5 


























0.5 



</STATIC_0PERAT0R_MATRIX> 
< /STATIC_0PERAT0R> 

######################################################## 

</SECTI0N-STATIC_0PERAT0RS> 

######################################################## 



######################################################## 
<SECTI0N-L0CAL_H0PPING_0PERAT0RS> 

######################################################## 

<L0CAL_H0PPING_0PERAT0R> 

# here comes the first local hopping operator 
<L0CAL_H0PPING_0PERAT0R_NAME> 

f _N+ 1 -dagger _up 



57 



< /LOCAL_HOPPING_OPERATOR_NAME> 
<LDCAL_HOPPING_OPERATOR_REPRESENTATION_INDEX> 
1 1 

</LOCAL_HOPPING_OPERATOR_REP RESENT ATION_INDEX> 
<LOCAL_HQPPING_OPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</LOCAL_HOPPING_OPERATOR_SIGN> 
<LOCAL_HOPPING_OPERATOR_MATRIX> 





10 

10 

</L0CAL_H0PPING_DPERAT0R_MATRIX> 
<LOCAL_DN_SITE_ENERGY_MATRIX> 





10 

1 

</LOCAL_DN_SITE_ENERGY_MATRIX> 
< /L0CAL_H0PPING_0PERAT0R> 

# here comes the second local hopping operator 
<L0CAL_H0PPING_0PERAT0R> 
<L0CAL_H0PPING_0PERAT0R_NAME> 

f _N+ 1 .dagger _down 
</L0CAL_H0PPING_0PERAT0R_NAME> 
<L0CAL_H0PPING_0PERAT0R_REPRESENTATI0N_INDEX> 
1 -1 

</L0CAL_H0PPING_0PERAT0R_REPRESENTATI0N_INDEX> 
<L0CAL_H0PPING_0PERAT0R_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</LQCAL_H0PPING_0PERATQR_SIGN> 
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<LOCAL_HOPPING_OPERATOR_MATRIX> 















1 





























-1 






</LOCAL_HOPPING_OPERATOR_MATRIX> 
<LOCAL_ON_SITE_ENERGY_MATRIX> 


















1 





























1 



< /LOCAL_ON_S ITE_ENERGY_MATRIX> 
< /LOCAL_HOPPING_OPERATOR> 

######################################################## 
</SECTION-LOCAL_HOPPING_OPERATQRS> 

######################################################## 
######################################################## 

<SECTION-LOCAL_HAMILTONIAN> 

# The total Local Hamiltonian of the added site 

# will be written as a sum of Local Hamiltonian terms 

# acting on site n: 

# The number of <LQCAL_HAMILTONIAN_TERM> items 

# must be equal with the variable ' local_coupling_no ' 

# set in the first section 

# The first tag is the name of the local Hamiltonian term 

# the second one gives the representation indices 

# last come the matrix elements . 

######################################################## 

# first term in the Local_Hamiltonian 

####################################################### 

<LOCAL_HAMILTONIAN_TERM> 

<LOCAL_HAMILTONIAN_NAME> 
H_l 

</LOCAL_HAMILTONIAN_NAME> 

<LOCAL_HAMILTONIAN_COUPLING> 

#irrelevant 
0.0 

</L0CAL_HAMILT0NIAN_C0UPLING> 
<L0CAL_HAMILT0NIAN_REP RESENT ATI0N_INDEX> 
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</LQCAL_HAMILTQNIAN_REPRESENTATIQN_INDEX> 
<LOCAL_HAMILTQNIAN_MATRIX> 



1 





























1 


















</LOCAL_HAMILTONIAN_MATRIX> 

</LOCAL_HAMILTONIAN_TERM> 

</SECTION-LOCAL_HAMILTONIAN> 

######################################################## 
<SECTIQN-SPECTRAL_FUNCTIQN> 

######################################################## 

# In this section we set up the spectral function 

# calculation for different types of operators. 

# All the time we provide in the input file the 

# matrices for the and B' operators 

# but the spectral function is computed for the 

# combination 

# spectral function for <f_0_up, f _0_dagger_up> 
{f _0_dagger_up; f_0 .dagger _up} 

# spectral function for <f_0_down, f _0_dagger_down> 
{f _0_dagger_down; f _0_dagger_down} 

#spectral function for <S_z, S_z> 
{S_z; S_z} 

######################################################## 

</SECTIQN-SPECTRAL_FUNCTION> 

######################################################## 

# The section below is not needed by the fnrg binary. 

# It is needed when doing the broadening of the spectral function 

# with the utility sfb. 

######################################################## 
<SECTION-SPECTRAL_FUNCTION_BROADENING> 

######################################################## 

# In this section we include broadening parameters that 

# are used for computing the broaden spectral function 

# The methods that can be used are either L0G_GAUSS or 

# INTERPOLATIVE_LOG_GAUSS . 
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# These procedures are described in the manual. 



# The quantum temperature is relevant only 

# for the INTERPOLATIVE_LOG_GAUSS method, 

# otherwise only the 'broadening_parameter ' parameter 

# setting the width of the log-normal 

# distribution function is relevant . 

# The quantum temperature must be always lower than T 

# for finite temperature calculations. 

# A good choice is to set the quatum temperature in the range 

# 0.5T - T. 

# The three parameters spectral_f unction_* 

# specify which spectral functions are being calculated. 

# The accepted values are YES or NO only. 

# The static_average_dmnrg flags specifies whether the expectation values 

# of the static operators are computed. 
#The accepted values are YES or NO only. 

# The green_function_dmnrg specifies whether the real part of the 

# Green's function is computed from the spectral function data. 

broadening_method = INTERPOLATIVE_LOG_GAUSS # method used either LOG.GAUSS or INTERPOLATIVE_LOG_GAUSS 
broadening_parameter = 0.70 # parameter used for the log gaussian distribution 
broadening_energy_minim = le-6 # minimum energy for which the broadening is performed 
broadening_energy_maxim = 4.0 # maximum energy for which the broadening is performed 
broadening_grid_mesh = 100 # grid mesh on a log scale between the minimum and maximum energies 
quantum.temperature = le-7 # used for the INTERP0LATIVE_L0G_GAUSS method only. 

spectral_function_nrg_even = YES # spectral function for the even iterations is computed 
spectral_function_nrg_odd = YES # spectral function for the odd iterations part is computed 
spectral_function_dmnrg = YES # spectral function for the dmnrg calculation is computed. 

static_average_dmnrg = YES # average functions similar with the case of spectral function. 

green_f unction_nrg_even = YES # calculation of the real/imaginary part of the Green's function 
from the even part of the nrg spectral function 

green_f unction_nrg_odd = YES # calculation of the real/imaginary part of the Green's function 
from the odd part of the nrg spectral function 

green_f unction_dmnrg = NO # calculation of the real/imaginary part of the Green's function 
from the spectral function 

green_f unction_grid_mesh = 400 # points for the mesh when doing Hilbert transform 

######################################################## 
</SECTI0N-SPECTRAL_FUNCTI0N_BR0ADENING> 

######################################################## 
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Appendix B 



Input file for the Kondo model 

with ^charge! 1 ) x SU spm (2) 

symmetries 



# The comments must start with the symbol # . 

# Empty lines are ignored. 

######################################################### 

<SECTION-PARAMETERS> 

######################################################### 

# Explanation of variables 

# 'model' represents the type of the model that is 

# used for the calculations. It can be any name. 

# 'lambda' represents the value for the logarithmic 

# discretization parameter. Usually it is fixed between 2 and 3. 

# 'max_state_no ' is the number of kept states 

# after an iteration. 

# ' iteration_no ' is the number of iterations performed 

# along the NRG run. Typically it is set to 40 - 70. 

# ' symmetry_no ' specifies how many symmetries 

# are used. The symmetries are characterized 

# in <SECTION-SYMMETRIES> 

# ' coupling_no ' specifies how many couplings 

# are present in the interaction Hamiltonian. 

# 'local_coupling_no' is the number of 

# Hamiltonians at the sites we add to the Wilson chain 

# 'spectral_operator_no' is the number of spectral operators 

# acting on site 0, we keep track of. 

# These spectral operators are characterized 
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# in section <SECTION-SPECTRAL-OPERATORS> 



# ' static_operator_no ' is the number of static operators 

# acting on site 0, we keep track of. 

# These static operators are characterized 

# in section <SECTI0N-STATIC-0PERAT0RS> 

# 'hopping_operator_no ' is the number of hopping operators. 

# These hopping operators are characterized 

# in section <SECTI0N-H0PPING-0PERAT0RS> 

# 'block_state_no ' specifies the number 

# of initial states/multiplets . These states are characterized 

# by the representation indices and are listed 

# in section <SECTI0N-BL0CK_STATES> 

######################################################### 
######################################################### 

model = kondojnodel # the name of the present model 
lambda = 2.0 # discretization parameter 
max_state_no = 100 # maximum number of kept multiplets 
iteration_no = 20 # allowed number of iterations 
symmetry_no = 2 # the number of symmetries 

coupling_no = 1 # no of Hamiltonian terms, i.e. J in this case 

spectral_operator_no = 2 # no of operators for which correlation functions are computed 
static_operator_no = # no of operators for which the static average is computed. 
hopping_operator_no = 1 # the number of hopping operators 
block_state_no = 4 # the number of initial block states 
local_state_no = 3 # the number of local states added in each iteration 
local_coupling_no = 1 # the number of Hamiltonian terms on the added site 
spectral_function_no = 2 # the number of spectral functions that need to be generated 
interval_no = 1000 # mash grid between the maximum and minimum energy at each iteration 

degeneracy_threshhold = le-6 # energy thresh-hold for discarding states 
temperature = 0.0 # temperature used for the run 

######################################################## 
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</SECTION-PARAMETERS> 

######################################################## 



######################################################## 

<SECTION-FLAGS> 

######################################################## 



# dmnrg_flag controls whether the backward procedure 

# has to be done and the full set of eigen-states be 

# used for the calculation of the spectral function. 

# ON - backwards iteration is done. 

# OFF - only the nrg calculation is performed. 

# text_swap_f iles_f lag = flag that fixes if the text mode is 

# used to save the files. The files will be saved in binary mode anyway. 

# ON - the files are saved in the text mode and kept on the disk. 

# 0FF= the files are not saved in the text mode. 

# binary _swap_f iles_f lag = flag that controls whether the 

# unnecessary binary files will be removed or not not after the calculation is done. 

# ON - the files will be kept on the disk. 

# OFF- the files will be removed from the disk. 

# hoppings_on_site_energies_f lag = flag that controls 

# the reading of the hoppings and the on-site energies. 

# ON - the hoppings and the on site energies 

# are read from the files in the results/mapping/ folder 

# These files are generated by the he utility. 

# OFF - the hoppings are computed on the fly 

# assuming a flat density of states (DOS =0.5) on [-1, 1], 

# and the on-site energies are set to zero. 

dmnrg_flag = ON 
text_swap_f iles_f lag = OFF 
binary_swap_f iles_f lag = ON 
hoppings_on_site_energies_f lag =0FF 



######################################################## 

< /SECTION-FLAGS > 

######################################################## 



######################################################### 

<SECTION-SYMMETRIES> 

######################################################### 

# As many symmetries must be defined below 

# as set by the variable ' symmetry_no ' 

# in the previous section. 

# The possible symmetry types are the following: 
#U(1), SU(2), Z(2), charge_SU(2) . 

# 
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# Next the limits for the representation index 

# are set . All states with smaller/larger 

# representation indices will be cut off 

# For U(l) symmetries the representations are labeled 

# by integers (2 S_z) and for SU(2) by non-negative 

# integers (2S) . 

# For the Z(2) symmetry the representations are labeled 

# by integers 

# - even parity states 

# 1 - odd parity states 

######################################################### 



# charge symmetry corresponding to Q 
U(l) -10 10 

# spin symmetry corresponding to S 
SU(2) 20 



######################################################### 

</SECTION-SYMMETRIES> 

######################################################### 



######################################################### 

<SECTI0N-BL0CK_STATES> 

######################################################### 

# This section gives the representation indices of the 

# block states of the first iteration. 

# The number of rows has to be equal 

# with the 'block_state_no ' in the first section, 

# and the number of columns to the ' symmetry_no ' . 

# Each column contains indices corresponding to a symmetry. 
######################################################### 



-1 


1 











2 


1 


1 



######################################################### 

</SECTION-BLQCK_STATES> 

######################################################### 



######################################################### 
< SECTI 0N-L0CAL_STATES > 

######################################################### 

# This section gives the representation indices 

# of the local states . The number of rows has to be equal 

# with the ' local_state_no ' in the first section 

# and the number of columns to the ' symmetry_no ' . 

# Each column contains indices corresponding to a symmetry. 
######################################################## 
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-1 

1 

1 

######################################################## 

</SECTION-LOCAL_STATES> 

######################################################## 



######################################################## 
<SECTION-LOCAL_STATES_SIGNS> 

######################################################## 

# These are the signs of the added states. 

# The number of signs has to be equal with 

# the number of ' local_state_no ' . 

# All the signs must be enumerated in one line . 

# For even electron number, the sign is 1, 

# while for odd electron number the sign is -1. 
######################################################### 

1-11 

######################################################### 
< /SECTION-LOCAL_STATES_SIGNS > 

######################################################### 



######################################################### 
<SECTION-BLOCK_HAMILTONIAN> 

# The total Hamiltonian will be written as a sum 

# of the Hamiltonian terms. 

rr \-^coupLirtg^no ^ rr 

Ho — l^s fastis 

# There has to be as many <HAMILTONIAN_TERM> items 

# as 'coupling_no' were set in the first section. 

# The first tag is the name of the Hamiltonian term, 

# the second one gives the representation indices and 

# at last comes the matrix elements that will be 

# read into the sectors. 



# first term in the Hamiltonian 

<BLOCK_HAMILTQNIAN_TERM> 

<BLOCK_HAMILTONIAN_NAME> 
HJCondo 

</BLOCK_HAMILTONIAN_NAME> 
<BLQCK_HAMILTQNIAN_CQUPLING> 

# J - exchange coupling 
J = 0.5 

< /BL0CK_HAMILT0NI AN_C0UPLING> 
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<BLOCK_HAMILTONIAN_REP RESENT ATION_INDEX> 


</BL0CK_HAMILTONIAN_REPRESENTATIQN_INDEX> 
<BLOCK_HAMILTONIAN_MATRIX> 


















-0.75 














0.25 


















</BLOCK_HAMILTONIAN_MATRIX> 
< /BLOCK_HAMILTONIAN_TERM> 

######################################################## 

</SECTION-BLOCK_HAMILTONIAN> 

######################################################## 
######################################################## 

<SECTI0N-H0PPING_0PERAT0RS> 

# There must be as many <H0PPING_0PERAT0R> items 

# as 'hopping_operator_no ' were set in the first section 

# The matrices must contain the Reduced Matrix Elements. 
######################################################## 

# here is the first hopping operator 
<H0PPING_0PERAT0R> 
<H0PPING_0PERAT0R_NAME> 

f _N -dagger 

</H0PPING_0PERAT0R_NAME> 
<H0PPING_0PERAT0R_REPRESENTATI0N_INDEX> 
1 1 

</H0PPING_0PERAT0R_REP RESENT ATI0N_INDEX> 
<H0PPING_0PERAT0R_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</H0PPING_0PERAT0R_SIGN> 
<H0PPING_0PERAT0R_MATRIX> 
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10 

10 

0.707106781186 -1.224744871392 

</HOPPING_OPERATOR_MATRIX> 



</H0PPING_0PERAT0R> 

######################################################## 
</SECTI0N-H0PPING_0PERAT0RS> 

######################################################## 



######################################################## 
<SECTION-SPECTRAL_OPERATORS> 

######################################################## 

# In this section we set up the matrices for the 

# operators for which the spectral function is computed 

# Matrices must contain the reduced matrix elements . 
####################################################### 

# here comes the first spectral operator. 



<SPECTRAL_0PERAT0R> 

<SPECTRAL_0PERAT0R_NAME> 

S 

</SPECTRAL_0PERAT0R_NAME> 

<SPECTRAL_0PERAT0R_REP RESENT ATI0N_INDEX> 
2 

</SPECTRAL_0PERAT0R_REPRESENTATI0N_INDEX> 
<SPECTRAL_0PERAT0R_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
1 

</SPECTRAL_0PERAT0R_SIGN> 
< SPECTRAL_0PERAT0R_MATRIX > 

0.866025403784 

-0.866025403784 

0.5 0.707106781186 



< / SPECTRAL_0PERAT0R_MATRIX> 
< /SPECTRAL_0PERAT0R> 

# here comes the second spectral operator. 
<SPECTRAL_0PERAT0R> 







0.866025403784 
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<SPECTRAL_OPERATOR_NAME> 



f_0 .dagger 

</SPECTRAL_OPERATOR_NAME> 

<SPECTRAL_OPERATOR_REP RESENT ATION_INDEX> 
1 1 

</SPECTRAL_OPERATOR_REPRESENTATIQN_INDEX> 
<SPECTRAL_OPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</SPECTRAL_OPERATOR_SIGN> 

< SPECTRAL_OPERATOR_MATRIX> 



1 
1 
0.707106781186 -1.224744871392 

< / SPECTRAL_0PERAT0R_MATRIX> 

< / SPECTRAL_0PERAT0R> 



######################################################## 
</SECTI0N-SPECTRAL_0PERAT0RS> 

######################################################## 



######################################################## 
<SECTI0N-L0CAL_H0PPING_0PERAT0RS> 

######################################################## 

<L0CAL_H0PPING_0PERAT0R> 

# here comes the first local hopping operator 
<L0CAL_H0PPING_0PERAT0R_NAME> 

f JJ+1 .dagger 

</L0CAL_H0PPING_0PERAT0R_NAME> 
<L0CAL_H0PPING_0PERAT0R_REPRESENTATI0N_INDEX> 
1 1 

</L0CAL_H0PPING_0PERAT0R_REP RESENT ATI0N_INDEX> 
<L0CAL_H0PPING_0PERAT0R_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
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-1 



</LOCAL_HOPPING_OPERATOR_SIGN> 
<LOCAL_HQPPING_QPERATOR_MATRIX> 



1 
-1.4142135623731 

</LOCAL_HOPPING_OPERATOR_MATRIX> 
<LOCAL_DN_SITE_ENERGY_MATRIX> 





</LOCAL_ON_SITE_ENERGY_MATRIX> 
< /LDCAL_HOPPING_OPERATOR> 



######################################################## 
</SECTION-LOCAL_HOPPING_OPERATORS> 

######################################################## 
######################################################## 

<SECTION-LOCAL_HAMILTONIAN> 

# The total Local Hamiltonian of the added site 

# will be written as a sum of Local Hamiltonian terms 

# acting on site n: 

tin — y \ C/i,nOi 7 n 

# The number of <L0CAL_HAMILT0NIAN_TERM> items 

# must be equal with the variable ' local_coupling_no ' 

# set in the first section 

# The first tag is the name of the local Hamiltonian term 

# the second one gives the representation indices 

# last come the matrix elements . 

######################################################## 

# first term in the Local_Hamiltonian 

####################################################### 

<L0CAL_HAMILT0NIAN_TERM> 

<LDCAL_HAMILTONIAN_NAME> 
H_l 

</L0CAL_HAMILT0NIAN_NAME> 

<L0CAL_HAMILT0NIAN_C0UPLING> 

#irrelevant 
0.0 
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< /LOCAL_HAMILTONI AN_COUPLING> 
<LOCAL_HAMILTONIAN_REP RESENT ATION_INDEX> 


</LOCAL_HAMILTONIAN_REPRESENTATIQN_INDEX> 
<LOCAL_HAMILTONIAN_MATRIX> 

1 

1 

</LOCAL_HAMILTONIAN_MATRIX> 

</LOCAL_HAMILTONIAN_TERM> 

</SECTION-LOCAL_HAMILTONIAN> 

######################################################## 
<SECTION-SPECTRAL_FUNCTION> 

######################################################## 

# In this section we set up the spectral function 

# calculation for different types of operators. 

# All the time we provide in the input file the 

# matrices for the and operators 

# but the spectral function is computed for the 

# combination 

# spectral function for <f_0, f_0_dagger> 
{f_0_dagger; f_0_dagger } 

#spectral function for spin <S, S> 
{S; S} 

######################################################## 

</SECTION-SPECTRAL_FUNCTION> 

######################################################## 



# The section below is not needed by the fnrg binary. 

# It is needed when doing the broadening of the spectral function 

# with the utility sfb. 

######################################################## 
<SECTION-SPECTRAL_FUNCTION_BROADENING> 

######################################################## 

# In this section we include broadening parameters that 

# are used for computing the broaden spectral function 

# The methods that can be used are either L0G_GAUSS or 

# INTERPOLATIVE_LOG_GAUSS . 

# These procedures are described in the manual. 
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# The quantum temperature is relevant only 

# for the INTERPOLATIVE_LOG_GAUSS method, 

# otherwise only the 'broadening_parameter ' parameter 

# setting the width of the log-normal 

# distribution function is relevant . 

# The quantum temperature must be always lower than T 

# for finite temperature calculations. 

# A good choice is to set the quatum temperature in the range 

# 0.5T - T. 

# The three parameters spectral_f unction_* 

# specify which spectral functions are being calculated. 

# The accepted values are YES or NO only. 

# The static_average_dmnrg flags specifies whether the expectation values 

# of the static operators are computed. 

# The accepted values are YES or NO only. 

# The green_function_dmnrg specifies whether the real part of the 

# Green's function is computed from the spectral function data. 

broadening_method = INTERPOLATIVE_LOG_GAUSS # method used either LOG.GAUSS or INTERPOLATIVE_LOG_GAUSS 
broadening_parameter = 0.70 # parameter used for the log gaussian distribution 
broadening_energy_minim = le-6 # minimum energy for which the broadening is performed 
broadening_energy_maxim = 4.0 # maximum energy for which the broadening is performed 
broadening_grid_mesh = 100 # grid mesh on a log scale between the minimum and maximum energies 
quantum.temperature = le-7 # used for the INTERP0LATIVE_L0G_GAUSS method only. 

spectral_function_nrg_even = YES # spectral function for the even iterations is computed 
spectral_function_nrg_odd = YES # spectral function for the odd iterations part is computed 
spectral_function_dmnrg = YES # spectral function for the dmnrg calculation is computed. 

static_average_dmnrg = YES # average functions similar with the case of spectral function. 

green_f unction_nrg_even = YES # calculation of the real/imaginary part of the Green's function 
from the even part of the nrg spectral function 

green_f unction_nrg_odd = YES # calculation of the real/imaginary part of the Green's function 
from the odd part of the nrg spectral function 

green_f unction_dmnrg = NO # calculation of the real/imaginary part of the Green's function 
from the spectral function 

green_f unction_grid_mesh = 400 # points for the mesh when doing Hilbert transform 

######################################################## 
</SECTI0N-SPECTRAL_FUNCTI0N_BR0ADENING> 

######################################################## 
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Appendix C 



Input file for the Kondo model 
with SU chdbTge (2) x SU spm (2) 
symmetries 

# The comments must start with the symbol # . 

# Empty lines are ignored. 

######################################################### 
< SECTI ON-P ARAMETERS > 

######################################################### 

# Explanation of variables 

# 'model' represents the type of the model that is 

# used for the calculations. It can be any name. 

# 'lambda' represents the value for the logarithmic 

# discretization parameter. Usually it is fixed between 2 and 3. 

# 'max_state_no ' is the number of kept states 

# after an iteration. 

# ' iteration_no ' is the number of iterations performed 

# along the NRG run. Typically it is set to 40 - 70. 

# ' symmetry_no ' specifies how many symmetries 

# are used. The symmetries are characterized 

# in <SECTION-SYMMETRIES> 

# ' coupling_no ' specifies how many couplings 

# are present in the interaction Hamiltonian. 

# 'local_coupling_no' is the number of 

# Hamiltonians at the sites we add to the Wilson chain 

# 'spectral_operator_no' is the number of spectral operators 

# acting on site 0, we keep track of. 

# These spectral operators are characterized 
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# in section <SECTION-SPECTRAL-OPERATORS> 



# ' static_operator_no ' is the number of static operators 

# acting on site 0, we keep track of. 

# These static operators are characterized 

# in section <SECTI0N-STATIC-0PERAT0RS> 

# 'hopping_operator_no ' is the number of hopping operators. 

# These hopping operators are characterized 

# in section <SECTI0N-H0PPING-0PERAT0RS> 

# 'block_state_no ' specifies the number 

# of initial states/multiplets . These states are characterized 

# by the representation indices and are listed 

# in section <SECTI0N-BL0CK_STATES> 

######################################################### 
######################################################### 

model = kondojnodel # the name of the present model 
lambda = 2.0 # discretization parameter 
max_state_no = 100 # maximum number of kept multiplets 
iteration_no = 20 # allowed number of iterations 
symmetry_no = 2 # the number of symmetries 

coupling_no = 1 # no of Hamiltonian terms, i.e. J in this case 

spectral_operator_no = 2 # no of operators for which correlation functions are computed 
static_operator_no = # no of operators for which the static average is computed. 
hopping_operator_no = 1 # the number of hopping operators 
block_state_no = 3 # the number of initial block states 
local_state_no = 2 # the number of local states added in each iteration 
local_coupling_no = 1 # the number of Hamiltonian terms on the added site 
spectral_function_no = 2 # the number of spectral functions that need to be generated 
interval_no = 1000 # mash grid between the maximum and minimum energy at each iteration 
degeneracy_threshhold = le-6 # energy thresh-hold for discarding states 
temperature = 0.0 # temperature used for the run 

######################################################## 
< / SECTI 0N-P ARAMETERS > 
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######################################################## 



######################################################## 

<SECTION-FLAGS> 

######################################################## 



# dmnrg_flag controls whether the backward procedure 

# has to be done and the full set of eigen-states be 

# used for the calculation of the spectral function. 

# ON - backwards iteration is done. 

# OFF - only the nrg calculation is performed. 

# text_swap_f iles_f lag = flag that fixes if the text mode is 

# used to save the files. The files will be saved in binary mode anyway. 

# ON - the files are saved in the text mode and kept on the disk. 

# 0FF= the files are not saved in the text mode. 

# binary _swap_f iles_f lag = flag that controls whether the 

# unnecessary binary files will be removed or not not after the calculation is done. 

# ON - the files will be kept on the disk. 

# OFF- the files will be removed from the disk. 

# hoppings_on_site_energies_f lag = flag that controls 

# the reading of the hoppings and the on-site energies. 

# ON - the hoppings and the on site energies 

# are read from the files in the results/mapping/ folder 

# These files are generated by the he utility. 

# OFF - the hoppings are computed on the fly 

# assuming a flat density of states (DOS =0.5) on [-1, 1], 

# and the on-site energies are set to zero. 

dmnrg_flag = ON 
text_swap_f iles_f lag = OFF 
binary_swap_f iles_f lag = ON 
hoppings_on_site_energies_f lag = OFF 

######################################################## 

< / SECTION-FLAGS > 

######################################################## 



######################################################### 

<SECTION-SYMMETRIES> 

######################################################### 

# As many symmetries must be defined below 

# as set by the variable ' symmetry_no ' 

# in the previous section. 

# The possible symmetry types are the following: 
#U(1), SU(2), Z(2), charge_SU(2) . 

# 

# Next the limits for the representation index 
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# are set . All states with smaller/larger 

# representation indices will be cut off 

# For U(l) symmetries the representations are labeled 

# by integers (2 S_z) and for SU(2) by non-negative 

# integers (2S) . 

# For the Z(2) symmetry the representations are labeled 

# by integers 

# - even parity states 

# 1 - odd parity states 

######################################################### 



# charge symmetry corresponding to Q 
charge_SU(2) 20 

# spin symmetry corresponding to S 
SU(2) 20 



######################################################### 

</SECTION-SYMMETRIES> 

######################################################### 



######################################################### 

<SECTI0N-BLOCK_STATES> 

######################################################### 

# This section gives the representation indices of the 

# block states of the first iteration. 

# The number of rows has to be equal 

# with the 'block_state_no ' in the first section, 

# and the number of columns to the ' symmetry_no ' . 

# Each column contains indices corresponding to a symmetry. 
######################################################### 



2 

1 1 

######################################################### 
</SECTION-BLQCK_STATES> 

######################################################### 



######################################################### 
< SECTI 0N-L0CAL_STATES > 

######################################################### 

# This section gives the representation indices 

# of the local states . The number of rows has to be equal 

# with the ' local_state_no ' in the first section 

# and the number of columns to the ' symmetry_no ' . 

# Each column contains indices corresponding to a symmetry. 
######################################################## 

1 

1 
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######################################################## 

</SECTION-LOCAL_STATES> 

######################################################## 



######################################################## 
<SECTION-LOCAL_STATES_SIGNS> 

######################################################## 

# These are the signs of the added states. 

# The number of signs has to be equal with 

# the number of ' local_state_no ' . 

# All the signs must be enumerated in one line . 

# For even electron number, the sign is 1, 

# while for odd electron number the sign is -1. 
######################################################### 

-1 1 

######################################################### 
< /SECTION-LOCAL_STATES_SIGNS > 

######################################################### 



######################################################### 
<SECTION-BLOCK_HAMILTONIAN> 

# The total Hamiltonian will be written as a sum 

# of the Hamiltonian terms. 

tj \-^couplinq^no /-i rr 

# There has to be as many <HAMILTONIAN_TERM> items 

# as 'coupling_no' were set in the first section. 

# The first tag is the name of the Hamiltonian term, 

# the second one gives the representation indices and 

# at last comes the matrix elements that will be 

# read into the sectors. 



# first term in the Hamiltonian 

<BLOCK_HAMILTONIAN_TERM> 

<BLOCK_HAMILTONIAN_NAME> 
HJCondo 

</BLOCK_HAMILTONIAN_NAME> 
<BLOCK_HAMILTONIAN_CQUPLING> 

# J - exchange coupling 
J = 0.5 

</BL0CK_HAMILT0NIAN_C0UPLING> 
<BL0CK_HAMILT0NIAN_REP RESENT ATI0N_INDEX> 
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</BLOCK_HAMILTONIAN_REPRESENTATION_INDEX> 



<BLQCK_HAMILTQNIAN_MATRIX> 



-0.75 






0.25 




</BLOCK_HAMILTONIAN_MATRIX> 
</BL0CK_HAMILT0NIAN_TERM> 

######################################################## 
< / SECTI0N-BL0CK_HAMILT0NI AN> 

######################################################## 
######################################################## 

<SECTIQN-HQPPING_OPERATORS> 

# There must be as many <H0PPING_0PERAT0R> items 

# as 'hopping_operator_no ' were set in the first section 

# The matrices must contain the Reduced Matrix Elements. 
######################################################## 

# here is the first hopping operator 
<H0PPING_0PERAT0R> 
<H0PPING_0PERAT0R_NAME> 

f _N_dagger 

</H0PPING_0PERAT0R_NAME> 
<HOPPING_QPERATOR_REPRESENTATION_INDEX> 
1 1 

</H0PPING_0PERAT0R_REP RESENT ATI0N_INDEX> 
<H0PPING_QPERAT0R_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

< /HOPP ING_0PERAT0R_S IGN> 
<H0PPING_0PERAT0R_MATRIX> 

-1.4142135623731 

-1.4142135623731 

0.707106781186 -1.224744871391 



SO 



</HOPPING_OPERATOR_MATRIX> 



</HOPPING_OPERATOR> 

######################################################## 
</SECTION-HOPPING_OPERATORS> 

######################################################## 



######################################################## 
<SECTION-SPECTRAL_OPERATORS> 

######################################################## 

# In this section we set up the matrices for the 

# operators for which the spectral function is computed 

# Matrices must contain the reduced matrix elements . 
####################################################### 

# here comes the first spectral operator. 

<SPECTRAL_OPERATOR> 

# here comes the spectral operator. 
<SPECTRAL_OPERATOR> 

<SPECTRAL_OPERATOR_NAME> 

f_0 .dagger 

</SPECTRAL_OPERATOR_NAME> 

<SPECTRAL_OPERATOR_REP RESENT ATION_INDEX> 
1 1 

</SPECTRAL_OPERATOR_REPRESENTATION_INDEX> 
<SPECTRAL_OPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</SPECTRAL_OPERATOR_SIGN> 
< SPECTRAL_OPERATOR_MATRIX> 






.707106781186 








1.224744871391 



-1.4142135623731 
-1.4142135623731 




</SPECTRAL_OPERATOR_MATRIX> 



< / SPECTRAL_0PERAT0R> 



# here comes the second spectral operator. 
<SPECTRAL_0PERAT0R> 
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<SPECTRAL_OPERATOR> 

<SPECTRAL_OPERATOR_NAME> 

S 

</SPECTRAL_OPERATOR_NAME> 

<SPECTRAL_OPERATOR_REP RESENT ATION_INDEX> 
2 

</SPECTRAL_OPERATOR_REPRESENTATIQN_INDEX> 
<SPECTRAL_OPERATOR_SIGN> 

# The sign is 1 for bosonic operators and -1 for fermionic operators. 
1 

</SPECTRAL_OPERATOR_SIGN> 

< SPECTRAL_OPERATOR_MATRIX> 

-0.866025403784 
0.5 0.707106781186 
0.866025403784 

</SPECTRAL_0PERAT0R_MATRIX> 

< / SPECTRAL_0PERAT0R> 



######################################################## 
</SECTI0N-SPECTRAL_0PERAT0RS> 

######################################################## 



######################################################## 
<SECTI0N-L0CAL_H0PPING_0PERAT0RS> 

######################################################## 
<L0CAL_H0PPING_0PERAT0R> 

# here comes the first local hopping operator 
<L0CAL_H0PPING_0PERAT0R_NAME> 

f JJ+l .dagger 

</L0CAL_H0PPING_0PERAT0R_NAME> 
<L0CAL_H0PPING_0PERAT0R_REPRESENTATI0N_INDEX> 
1 1 

</L0CAL_H0PPING_0PERAT0R_REPRESENTATI0N_INDEX> 
<L0CAL_H0PPING_0PERAT0R_SIGN> 
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# The sign is 1 for bosonic operators and -1 for fermionic operators. 
-1 

</LOCAL_HOPPING_OPERATOR_SIGN> 
<LOCAL_HQPPING_QPERATOR_MATRIX> 

1.4142135623731 
-1.4142135623731 

</LOCAL_HOPPING_OPERATOR_MATRIX> 
<LOCAL_DN_SITE_ENERGY_MATRIX> 




</LOCAL_ON_SITE_ENERGY_MATRIX> 
< /L0CAL_H0PPING_0PERAT0R> 



######################################################## 
</SECTION-LOCAL_HDPPING_OPERATDRS> 

######################################################## 
######################################################## 
<SECTION-LOCAL_HAMILTONIAN> 

# The total Local Hamiltonian of the added site 

# will be written as a sum of Local Hamiltonian terms 

# acting on site n: 

tin — y \ C/i,nOi 7 n 

# The number of <LQCAL_HAMILTONIAN_TERM> items 

# must be equal with the variable ' local_coupling_no ' 

# set in the first section 

# The first tag is the name of the local Hamiltonian term 

# the second one gives the representation indices 

# last come the matrix elements . 

######################################################## 

# first term in the Local_Hamiltonian 

####################################################### 

<LOCAL_HAMILTONIAN_TERM> 

<L0CAL_HAMILTONIAN_NAME> 
H_l 

</LOCAL_HAMILTONIAN_NAME> 

<LOCAL_HAMILTONIAN_COUPLING> 

#irrelevant 
0.0 

</L0CAL_HAMILT0NIAN_C0UPLING> 
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<LOCAL_HAMILTONIAN_REP RESENT ATION_INDEX> 


</LQCAL_HAMILTQNIAN_REPRESENTATION_INDEX> 
<LOCAL_HAMILTONIAN_MATRIX> 

1 


</LOCAL_HAMILTONIAN_MATRIX> 
</LOCAL_HAMILTONIAN_TERM> 
< / SECTION-LOCAL_HAMILTONI AN> 

######################################################## 
<SECTION-SPECTRAL_FUNCTION> 

######################################################## 

# In this section we set up the spectral function 

# calculation for different types of operators. 

# All the time we provide in the input file the 

# matrices for the A* and operators 

# but the spectral function is computed for the 

# combination <A,B*> 

# spectral function for <f_0, f_0_dagger> 
{f_0_dagger; f_0_dagger } 

# spectral function for spin <S, S> 
{S; S } 

######################################################## 

</SECTION-SPECTRAL_FUNCTION> 

######################################################## 

# The section below is not needed by the fnrg binary. 

# It is needed when doing the broadening of the spectral function 

# with the utility sfb. 

######################################################## 
<SECTION-SPECTRAL_FUNCTION_BROADENING> 

######################################################## 

# In this section we include broadening parameters that 

# are used for computing the broaden spectral function 

# The methods that can be used are either L0G_GAUSS or 

# INTERPOLATIVE_LOG_GAUSS . 

# These procedures are described in the manual. 

# The quantum temperature is relevant only 

# for the INTERPOLATIVE_LOG_GAUSS method, 
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# otherwise only the 'broadening_parameter ' parameter 

# setting the width of the log-normal 

# distribution function is relevant . 

# The quantum temperature must be always lower than T 

# for finite temperature calculations. 

# A good choice is to set the quatum temperature in the range 

# 0.5T - T. 

# The three parameters spectral_f unction_* 

# specify which spectral functions are being calculated. 

# The accepted values are YES or NO only. 

# The static_average_dmnrg flags specifies whether the expectation values 

# of the static operators are computed. 

# The accepted values are YES or NO only. 

# The green_function_dmnrg specifies whether the real part of the 

# Green's function is computed from the spectral function data. 

broadening_method = INTERPOLATIVE_LOG_GAUSS # method used either LOG.GAUSS or INTERPOLATIVE_LOG_GAUSS 
broadening_parameter = 0.70 # parameter used for the log gaussian distribution 
broadening_energy_minim = le-6 # minimum energy for which the broadening is performed 
broadening_energy_maxim = 4.0 # maximum energy for which the broadening is performed 
broadening_grid_mesh = 100 # grid mesh on a log scale between the minimum and maximum energies 
quantum.temperature = le-7 # used for the INTERP0LATIVE_L0G_GAUSS method only. 

spectral_function_nrg_even = YES # spectral function for the even iterations is computed 
spectral_function_nrg_odd = YES # spectral function for the odd iterations part is computed 
spectral_function_dmnrg = YES # spectral function for the dmnrg calculation is computed. 

static_average_dmnrg = YES # average functions similar with the case of spectral function. 

green_f unction_nrg_even = YES # calculation of the real/imaginary part of the Green's function 
from the even part of the nrg spectral function 

green_f unction_nrg_odd = YES # calculation of the real/imaginary part of the Green's function 
from the odd part of the nrg spectral function 

green_f unction_dmnrg = NO # calculation of the real/imaginary part of the Green's function 
from the spectral function 

green_f unction_grid_mesh = 400 # points for the mesh when doing Hilbert transform 

######################################################## 
</SECTI0N-SPECTRAL_FUNCTI0N_BR0ADENING> 

######################################################## 
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Appendix D 

License agreements 



GNU LESSER GENERAL PUBLIC LICENSE Version 3, 29 June 2007 

Copyright (C) 2007 Free Software Foundation, Inc. <http://fsf.org/> Everyone is permitted 
to copy and distribute verbatim copies of this license document, but changing it is not allowed. 
This version of the GNU Lesser General Public License incorporates the terms and conditions 
of version 3 of the GNU General Public License, supplemented by the additional permissions 
listed below. 

0. Additional Definitions. 

As used herein, "this License" refers to version 3 of the GNU Lesser General Public License, 
and the " GNU GPL" refers to version 3 of the GNU General Public License. 
"The Library" refers to a covered work governed by this License, other than an Application or 
a Combined Work as defined below. 

An " Application" is any work that makes use of an interface provided by the Library, but which 
is not otherwise based on the Library. Defining a subclass of a class defined by the Library is 
deemed a mode of using an interface provided by the Library. 

A "Combined Work" is a work produced by combining or linking an Application with the 
Library. The particular version of the Library with which the Combined Work was made is 
also called the " Linked Version" . 

The "Minimal Corresponding Source" for a Combined Work means the Corresponding Source 

for the Combined Work, excluding any source code for portions of the Combined Work that, 

considered in isolation, are based on the Application, and not on the Linked Version. 

The "Corresponding Application Code" for a Combined Work means the object code and/or 

source code for the Application, including any data and utility programs needed for reproducing 

the Combined Work from the Application, but excluding the System Libraries of the Combined 

Work. 

1. Exception to Section 3 of the GNU GPL. 

You may convey a covered work under sections 3 and 4 of this License without being bound by 
section 3 of the GNU GPL. 

2. Conveying Modified Versions. 

If you modify a copy of the Library, and, in your modifications, a facility refers to a function or 
data to be supplied by an Application that uses the facility (other than as an argument passed 
when the facility is invoked), then you may convey a copy of the modified version: 

a) under this License, provided that you make a good faith effort to ensure that, in the event 
an Application does not supply the function or data, the facility still operates, and performs 
whatever part of its purpose remains meaningful, or 

b) under the GNU GPL, with none of the additional permissions of this License applicable to 
that copy. 

3. Object Code Incorporating Material from Library Header Files. 
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The object code form of an Application may incorporate material from a header file that is part 
of the Library. You may convey such object code under terms of your choice, provided that, 
if the incorporated material is not limited to numerical parameters, data structure layouts and 
accessors, or small macros, inline functions and templates (ten or fewer lines in length), you do 
both of the following: 

a) Give prominent notice with each copy of the object code that the Library is used in it and 
that the Library and its use are covered by this License. 

b) Accompany the object code with a copy of the GNU GPL and this license document. 

4. Combined Works. 

You may convey a Combined Work under terms of your choice that, taken together, effectively 
do not restrict modification of the portions of the Library contained in the Combined Work 
and reverse engineering for debugging such modifications, if you also do each of the following: 

a) Give prominent notice with each copy of the Combined Work that the Library is used in it 
and that the Library and its use are covered by this License. 

b) Accompany the Combined Work with a copy of the GNU GPL and this license document. 

c) For a Combined Work that displays copyright notices during execution, include the copyright 
notice for the Library among these notices, as well as a reference directing the user to the copies 
of the GNU GPL and this license document. 

d) Do one of the following: 

0) Convey the Minimal Corresponding Source under the terms of this License, and the Corre- 
sponding Application Code in a form suitable for, and under terms that permit, the user to 
recombine or relink the Application with a modified version of the Linked Version to produce a 
modified Combined Work, in the manner specified by section 6 of the GNU GPL for conveying 
Corresponding Source. 

1) Use a suitable shared library mechanism for linking with the Library. A suitable mechanism 
is one that (a) uses at run time a copy of the Library already present on the user's computer 
system, and (b) will operate properly with a modified version of the Library that is interface- 
compatible with the Linked Version. 

e) Provide Installation Information, but only if you would otherwise be required to provide such 
information under section 6 of the GNU GPL, and only to the extent that such information 
is necessary to install and execute a modified version of the Combined Work produced by 
recombining or relinking the Application with a modified version of the Linked Version. (If 
you use option 4d0, the Installation Information must accompany the Minimal Corresponding 
Source and Corresponding Application Code. If you use option 4dl, you must provide the 
Installation Information in the manner specified by section 6 of the GNU GPL for conveying 
Corresponding Source.) 

5. Combined Libraries. 

You may place library facilities that are a work based on the Library side by side in a single 
library together with other library facilities that are not Applications and are not covered by 
this License, and convey such a combined library under terms of your choice, if you do both of 
the following: 

a) Accompany the combined library with a copy of the same work based on the Library, un- 
combincd with any other library facilities, conveyed under the terms of this License. 

b) Give prominent notice with the combined library that part of it is a work based on the 
Library, and explaining where to find the accompanying uncombincd form of the same work. 

6. Revised Versions of the GNU Lesser General Public License. 

The Free Software Foundation may publish revised and/or new versions of the GNU Lesser 
General Public License from time to time. Such new versions will be similar in spirit to the 
present version, but may differ in detail to address new problems or concerns. 
Each version is given a distinguishing version number. If the Library as you received it specifies 
that a certain numbered version of the GNU Lesser General Public License " or any later version" 
applies to it, you have the option of following the terms and conditions either of that published 
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version or of any later version published by the Free Software Foundation. If the Library as 
you received it does not specify a version number of the GNU Lesser General Public License, 
you may choose any version of the GNU Lesser General Public License ever published by the 
Free Software Foundation. 

If the Library as you received it specifics that a proxy can decide whether future versions of the 
GNU Lesser General Public License shall apply, that proxy's public statement of acceptance of 
any version is permanent authorization for you to choose that version for the Library. 
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